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NACA  REPORT  12^5 

By  Morris  Morduchow 
1955 

Pa^e  l4,  column  2,  paragraph  2  of  section  entitled  "Separation,"  line  6 
should  read: 

"both  ^sep  (eq.  (65))  and  Ag(|)  (eq.  (64)).  However,  numerical" 


Page  19,  table  II:  The  values  2,590  and  38I  in  the  bottom  row  shoiild 

be  realined  to  indicate  that  the  multiplying  factor  "x  10-^"  in  the 
first  row  applies  also  to  them.  Thus,  the  table  should  appear  as 
follows ; 


Flow 

Rco,  cr  for 

I-loo  of  — 

0 

1 

^1/^00  =  1-5 

350  X  lo5 

104  X  10^ 

(adverse  pressure  gradient;  data  from  ref.  2) 

^l/u«  =  1 

2,390 

581 

(zero  pressiure  gradient;  data  from  ref.  1) 

Page  19,  table  III:  The  minus  sign  in  the  subtitle  of  part  (b)  shovuLd 
be  an  equal  sign,  so  that  this  subtitle  should  read: 

(b)  Values  over  flat  plate;  =  |  =0.8 
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Ky  Mokims  Mohim  <‘mow 


SUMMARY 

.1  ■'Oirrti/  Ilf  iiileijnil  iiaf)iiiih  in  litniiniir-hiinnilni ii-lann- 
itniihjsis  is  first  iiirin.  ,1  sniijtlr  inn/  siijficientlij  iirciiritlr 
nii  thiiil  fur  jirilcticiil  jiiiriiosi  s  of  l•<ll(•llllltlnl|  thi  jniifit  rfn  s 
(inclllilinil  sliibilifll\  Ilf  tin  lininnilr  cnlli  jni  ssibli  bninnlilrii 
lili/ir  in  iin  ii.riiil  jirissiirr  ijntiliint  irilb  hull  trnnsfrr  at  t/n- 
iralt  is  tin  n  jiiTsi  ntnl.  Fur  jhnr  nrir  a  flat  /ilatt.  the  nn  tlioil 
is  a jijilieablf  for  an  arbifrarili/  /iresrribeil  ihstribation  nf 
lein jieratarf  alonii  tin  siirfave  ami  for  ani/  ijiren  cansfant 
Franilll  narnber  c/nse  In  inii/>/.  Fur  flmr  m  a  /iressare  ijraiinnf. 
tin  nn-fhiiil  is  basfil  mi  a  Pramltl  nuinbir  nf  iinili/  'inil  a  anifnrm 
mail  teinperalnre.  simple  ami  acvarate  mellnnl  nf  ilelermin- 
iinj  the  separatinn  pnint  in  a  cnmpressible  Jlmr  mith  an  ailrerse 
pressure  ijrailient  nrir  a  surface  at  a  ijiren  unifnnn  mall  temper¬ 
ature  is  ilereinjieil .  I  he  anah/sis  is  haseil  nn  an  extensmn  nf 
the  Kdrnnin-Fnhihausen  melhml  tn  the  mnmentum  am!  thermal 
inerijii  equal  inns  in  cnnjunctinn  mith  fnurth-  ami  esperiallp 
Imther  ileijiee  rehmitii  ami  stintnatinn-e ntha! pq  jirnjiles.  Fmm 
the  equatinns  ihrirnl  hire,  cnnii usimis  reitiirihnii  the  effect  nf 
pressure  iiraiHent,  Mach  number,  ami  mall  temperature  nn  the 
hnumlarq-layer  characteristics  are  ilerireil  ami  i/lustrateil.  In 
particular  the  e  ffects  nn  ski n-frictinn .  heat-transfer  cnefficient. 
SI  jiaratinn  pnint  in  an  ailrerse  pressure  iirailient.  ami  stahililq 
nf  the  laminar  hnumlarq  lai/i  r  are  anali/:eil. 

INTKODK  TiON 

The  [)uri)()sr  of  tlic  ()r(‘S('nt  i'<‘[)ort  is  to  prosont  u  coiuprc- 
hciisivf  simimary  of  tiicoictical  invest  i>ra(  ions  of  eoiiipres- 
sible  laminar  bonndai  v  layers  wliieli  liavi'  been  <-arric'<l  out 
since  1940  at  the  I’olyteehnie  Institute  of  Ibooklyn  iimler 
tlie  sponsorslii])  and  witli  tlie  financial  assistance  of  the 
National  Advisory  Committee  for  Aeronanties.  The  results 
of  these  invn'stioations  are  <’ontained  primarily  in  refenaiees 
1  to  7.‘ 

briefly,  reference  1  is  an  inves(i<i;alion  of  the  rela(iv<>  iiierils 
of  various  types  of  intcfiral  methods  for  the  analysis  of 
laminar  boundary  la\-ers.  It  is  eonehided  that  the  one- 
|)arameler  method  base<l  on  the  Karmiin  momentum  integral 
equation  in  eonjntietion  with  sixih-dc'frree  velocity  profiles 
api)ears.  on  the  whole,  to  be  the  most  jiromisinf;  method  for 
analyziiifr  laminar  boundary  layers  in  oein'ral.  On  the  basis 
(»f  this  conclusion,  a  .simple  method  of  calculatin}'  compre.s- 
sible-bonndary-la.\ er  characteristics,  inclttdino  .separation 
f)oint  aird  stability  characteristics,  in  flow  with  a  pre.ssiire 
firadient.  witliotil  heat  tianslei'  at  the  wall  is  d(‘velopeil  in 


I  reb'rence  2.  It  is  further  shown  in  reference  2  that  iourth- 
i  <b‘trree  pi'oliles  are  ])referable  foi'  analyzing  si  agnation  Hows, 
ahile  the  se/)aialion  jjoini  in  an  adverse  pjessure  pradieiil 
I  I'an  be  still  more  acenratel\'  picdieled  with  se\ cm  h(lep:i  ee 
I  veloeitv  proliles.  In  refel'cnce  2.  it  is  further  verifieil,  by 
I  considering  the  flow  o\cr  a  Ihit  plate  with  heal  liansfer. 

I  that  sixt h-defxree  profih's  yield  lesnlts  of  sullii'iciil  aeeunie\ 

I  for  stability  ealeulat  ions.  In  I  he  eouise  of  such  ealeulat  ions, 
j  certain  modifications  in  the  approximate  stability  eriteiia  <if 
reference  S  wei'e  made,  and  these  are  shown  briell\-  in  refer¬ 
ence's  to  whii'h  are  essi'iit iall.\  a  summain  of  an  unpub- 
I  lished  report  by  I’rofessor  M.  Bloom  of  thi'  Bolyti'clinie 
'■  Institute  of  Brool\l\  ii  entitled  "( 'alenlal ion  of  Stabilit  e  (>f 
I  ( 'onstant-l'ressnre  Boundaiy  La.\crson  Isothermal  Surfaces 
j  With  an  Inteoral-Methoil  Mean-Flow  Solution.’'  This  report 
I  is  availitble  for  loan  or  reference  in  the  Division  of  Keseareh 
i  Information.  National  .\dvisor.\  ('ominittee  for  .Aeronaut ies, 

!  Washinptton.  I).  ('.  In  refereni’e  l>.  a  method  of  ealeiilatint^ 
1  the  comprt'ssible  laminar  bomidar.v  layer  in  a  pressure  oya- 
I  dient  w  ith  heat  t ransfer  is  developi'd .  This  reference  ineludi's 
till'  ealeilliilion  of  the  houndar.e  la.ver  over  a  flat  [)late  (zero 
;  pressure  frrailient )  w  ith  a  nonntiiform  wall  temperattux'.  Kef- 
I  ereiiee  7.  finally,  a|>plies  the  ('(plat ions  (leV('lop('(l  in  r('f('ren(  (' 

I  ti  to  a  freneral  study  of  tin'  ('fr('et  of  prt'ssuix'  gradient,  wall 
I  t('mp('rat nix',  and  Mach  nundx'r  on  t In' skin-friet ion,  lu'iit- 
!  transb'i'.  s('|)arat ion,  and  stability  eharaet('risties  of  laminar 
i  bonndaiw  lay('rs.  A  nu'thod  of  ealeulat in>r  tlu'  s('paration 
j  point  in  an  iiilvcrsi'  pix'ssure  ^rradii'iil  with  heal  transfer  is 
inclinh'd  tlu'ix'.  .Nimu'rieal  ('.xamiih's  ar('  also  iiiehuh'd  in 
r('f('r('ne('  7  to  ilhlslrat('  in  (h'tinl  tin'  eonelnsions  |■('a(■h('(l 
tlu'iX'.  riu'  nu'lhods  of  r('f('r('ne('s  ti  and  7  art'  ('Xt('nsions  of 
tin'  corn'spondino  nu'thods  of  r('f('r('ne('  2  to  eas('s  of  lu'al 
transb'r  at  the  wall. 

Tilt'  ('injihasis  in  tin'  pr('S('nt  re])ort  will  Ix'  on  tlu'  dev('lop- 
iiK'iit  of  nn'thods  of  ealcidation  of  laminar-compix'ssibh'- 
boundar.v-layt'r  characteristics.  In  particular,  a  im'thod  of 
calculatin';  the  boundary  layer  over  a  flat  plat('  with  a  non- 
isotlu'rmal  stirfaei',  that  is,  with  a  piivt'n  distribution  of 
ti'inperatinx'  along  the  wall,  will  be  pri'senti'd.  Tin'  IVandtl 
iiumln'r,  although  constant ,  is  h'ft  arbitrary  but  must  bt'ofthe 
orderof  magnitude  of  unit.v.  For  flow  in  a  pr('ssur(' gra(li('nt , 
a  iiK'thod  isgiv('n  of  calculating  tin'  boundary  layer  for  a  given 
distribution  of  xu'hx'ity  outside  the  boundary  layt'r,  a  gixx'ii 
ri'fi'ri'in'c  Nfat'h  ntimber,  and  a  given  iinirorm  wall  teinpi'ia- 
i  turi'.  'J'he  I’randtl  number  is  now  assumed  to  Ix'  tmity.  A 


'  Simt'  •>!  iln'sc  rcfiTruct's  ntnlnifis  tipjirfcitUilo  of  hitorost  Hot  in  th«‘  prosoni  n*|M»r!.  t!i»*  latlor  dots  not  aciijiilly  siijuTsotlr  any  of  those  roferoiKT.s. 
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liiclliod  of  calciilal iiifr  the  si'imial ion  point  in  coiiipiossihlo 
flow  with  an  advcisc  iiicssnic  gradient  over  a  surface  at  a 
given  iiniforni  wall  teinjieratnre  is  also  presented.  'I'lie 
niatlieiiiatieal  analysis  on  which  tlu‘se  methods  arc  based  will 
he  given  in  sufficient  rletail  to  show  elearly  the  logical  diweloj)- 
ment  of  tiu'  methoils  and  the  main  a]>pro.vimaling  assump¬ 
tions  which  are  made  as  well  as  the  range  of  applicability  of 
the  methods. 

In  addition  to  the  methoils  of  ealeulation,  the  implications 
of  the  equations  developi'd  hiu'e  regarding  tlie  various  char- 
aeteristies  of  laminar  boundary  layers  as  they  are  affecteii 
by  paraiueti'is  such  as  the  pressure  grailient,  wall  tempera¬ 
ture,  and  Mach  number  will  be  iliscussed.  The  discu.ssion 
will  include  stability  characteristics.  A  variety  of  mimiu'ical 
e.xamples  will  also  be  iliseussed.  and  details  of  a  few  of  these 
will  be  presented  here.  Further  di'tails  of  these  examples 
can  be  found  in  referenees  1  to  7. 

It  may  be  remarked  that,  as  is  well  known,  the  lili'ratiire 
on  laminar  boutidarv  layers  and  related  |)roblems  is  exten¬ 
sive  and  rich.  Ituleed,  wherever  pertinent,  reference  will 
be  included  in  this  report  to  ri'cent  work  which  has  appeared 
either  more  or  less  simultatieotisly  with,  or  since,  references 
1  to  7.  The  advantage  of  the  methods  of  calculation  devel- 
oped  in  this  report  is  that  they  combine  the  merits  of  ade¬ 
quate  accuracy  and  relative  eas(>  of  calculation.  (An  ordi¬ 
nary  desk  calculator  will  be  found  to  be  more  thatt  adequate 
for  all  of  the  ealculatiotis.  In  fact,  a  large  nttmber  of  the 
calculations  may  even  be  performed  by  a  stamiard  slide 
rule.  I  The  mathematical  analysis  will  likewise  entail  prac¬ 
tically'  the  minitiium  number  of  approximatittg  assumptions 
required  to  retain  both  simjilicity  atid  adi'qitate  accuracy. 

With  res|)('ct  to  the  pertinent  literature,  it  will  suflici',  at 
this  jioint,  to  tnention  brielly  theoretical  invi'stigatiotis  on 
the  general  case  of  the  compressible  laminar  boundary  layer 
with  pressure  gradient  and  heat  transfer.  Oidy  a  very 
limitod  number  of  exact  solutions,  that  is,  solutions  (which 
may  be  numerical)  based  on  solving  directly  the  original 
partial  ililTerential  eipiations  of  the  boundary  layer  es,sen- 
tially  without  any  mathematical  apiuoximations,  n|>pear  to 
have  been  obtained  thus  far.  Such  exact  solutions  are 
restricted  to  particular  types  of  (lows.  For  low-speed  (zero 
Mach  number),  but  tievertheless  compressibh*,  flows  with 
heat  transfer,  numerical  solutions  for  “wedge  flows,’’  or 
(lows  in  which  the  velocity  outside  the  boundary  layer  is 
pro|)ortional  to  a  power  of  the  axial  distance,  have  been 
developed  in  references  !l  to  11.  These  solutions  include  a 
small  normal  ma.ss  flow  at  the  wall.  Such  solutions  have 
been  recently  extended  for  nonuegligible  Mach  numbers  to 
cases  in  which  the  local  Mach  number  outside  of  the  bound¬ 
ary  layei'  is  pro|)ortional  to  a  (lower  of  the  axial  distance  (ref. 
12).  \  class  of  similar  solutions  (refs.  13  to  15)  for  high- 

s()eed  (i.  e.,  nonzero  Mach  number)  flows  has  been  rei'entlv 
derived,  under  certain  conditions,  and  calculated  with  the 
aid  of  electrotiic  computors. 

In  addition  to  referimces  ti  and  7,  several  a])i)roxiinale 
analyses  of  laminar  eomitressible  boundary  layers  with  pres¬ 
sure  gradient  and  heat  transfer  have  been  made.  An  analy¬ 
sis.  for  example,  based  on  a  ty|)e  of  apitroximatioti  used  bv 
Lighthill  has  heen  recently  made  in  referenci*  10,  with  empha¬ 


sis  on  se|)aiation.  \  method  based  on  "internal  "  and 
“externar'  solutions  of  the  compressible-boundary-layer 
eijuations,  (Jieviously  introduced  by  Karinan  and  Millikan 
in  an  analysis  for  incom()iessible  (low,  is  develo|)ed  in  refi'i- 
ence  17.  Keference  l.S  (Uesentsa  method  based  on  solutions 
which  have  beim  obtained  for  wedge  Hows.  References  l!i 
to  21  <leveh)|)  methods  based  on  an  extension  of  the  Karnnin- 
Fohlhausen  method  (with  fouith-degree  velocity  and  stag- 
nation-enthalin  (iroliles)  to  the  theiiual-energy,  as  well  as 
the  monuaitum,  (tartial  did'erential  eipialion.  (Kefs.  is.  2li. 
and  21  include  the  case  of  a  small  normal  mass  (low  at  the 
wall.)  Further  analyses  based  on  integral  Jiiethods  are 
given  in  references  22  to  24.  the  latter  being  a  study  of  the 
heat-insulating  |)ro|)erl.ies  of  the  laminar  boundary  layer,  .V 
small-|>erturbation  ty|)e  of  analy  sis  is  deve.lo|Je(l  in  reference 
25.  It  may  be  noted  that  all  of  the  foregoing  a|)|)roximaie 
analyses,  with  the  exce|)tion  of  references  Iti  to  IS  anil  25. 
are  ba.sed  on  integral  methods. 

'I'he  analysis  in  the  (treseiit  re|)ort  is,  for  sim|)licily,  based 
on  the  usual  assum|)tion  of  constant  s|)ecific  heats.  A  means 
of  taking  into  account  variable  s|)eci(ic  heat  with  tempera¬ 
ture,  at  least  for  (low  over  a  Hal  (ilate,  is  disenssed,  lor 
examide,  in  referenees  2t>  and  27.  As  has  already  been 
stall'd,  it  is  further  assumed  in  this  reitort  that  for  How  w  iili 
a  (u'cssure  gradient  the  I’randll  number  is  1.  .\n  a|)i)io\i- 

mate  means,  for  the  ease  of  zero  heat  transfer,  of  taking  into 
account  a  Prandtl  number  iliH'erent  from  unity  is  discussed  in 
reference  28  and  a|>i)lied  in  reference  2!t.  In  the  case  of  heat 
transfer  at  the  wall,  an  ai)|>roximate  means  of  taking  into 
account  a  Prandtl  number  other  than  1  would  be  to  multi|)ly 
the  Nusselt  number  (i.  e.,  heat-transfer  eoeflieient  J  obtained 
in  accordanee  with  the  method  given  here  by  the  cube  root 
of  the  Prandtl  number  (ef..  e.  g.,  refs.  30  and  31  i.  Siieti  a 
correction,  however,  may  be  considerably  inaccurate  at  very 
high  Mach  numbers  (ref,  12).  A  further  assumiition  in  the 
(jresent  analysis  of  How  irlth  a  jinxxiin  iinlilii  ut  is  that  the 
wall  tenqierature  is  uniform.-  A  summary  of  investigations 
on  How  over  a  nonisothermal  surface  in  a  (iressure  gradient 
(as  well  as  over  a  Hat  |)late)  is  given  in  reference  32.  Further 
information  can  also  be  obtained  in  reference  33.  Finally, 
it  must  be  noted  that  the  (tresent  investigation  is  based 
on  the  assum|)tion  that  the  coelfieient  of  viscosity  is  |)ro- 
|)ortional  to  the  absolute  tem|)erature,  with  the  (troportion- 
ality  factor  determined  so  that  Sutherland’s  relation  is 
exactly  satisfied  at  the  wall.  ’I'his  is  an  assumittion  com¬ 
monly  made  (cf.  ref.  34)  to  sim|)lify  the  analysis  and  yet 
retain  the  main  actual  influence  of  the  de])endence  of  the 
viscosity-  eoeflieient  on  temiierature,  at  least  for  .Mach 
numbers  below  5. 

I’he  present  re[)ort  is  divided  into  live  main  sections,  'riie 
first  section  discusses  concisely  the  various  main  types  of 
integral  mi'thods  in  laminar-hoimdary-layer  analysis  and 
their  relative  merits,  'I'he  second  section  develo|)s  the 
basic  equations  to  be  used  in  the  present  analysis.  'Phese 
eqtiations  are  valid  for  an  arbitrary  constant  Prandtl  tuimber 
(close  to  unity)  and  a  nonutiiform  wall  temperature.  In  the 


*  I'lip  pquiHliins,  howpvpr,  tlovclttpcl  in  llu‘  srotitiii  “  Husic  K<i nations,”  wliicli 

lain  lo  flows  wiilt  or  wiihoiil  a  prpsstire  urailiont.  arp  vali'l  als<'  for  noniiniforni  wall 
tPinjH'ratun*. 
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tliiril  section,  tliese  (‘quations  ai'c  ajiplied  (o  |>i'i‘sent  a  ■ 
nielliod  for  llie  culculalioii  of  the  t)onndary  layer  over  a 
Hat  plate  with  a  preserihed  distrihutioii  of  widl  temperature 
while  the  Prandti  number  is  kept  arbitrary.  In  tin'  fourth 
.section,  the  basic  eiiuations  are  used  to  yield  a  method  of 
ealenlatinn  the  boundary  layer  in  a  i;iven  jiressure  gradient 
over  a  surface  at  a  prescribed  uniform  wall  t('mperature. 
Here  the  Prandti  number  is  assumed  as  unity.  'I'he  ealeula- 
tion  of  the  separation  point  in  an  adverse  prt'ssure  Gradient 
is  inchuled  in  this  section.  'I'he  fifth  section,  finally,  dis¬ 
cusses  the  various  f;en('ral  conclusions  on  the  boundary -layer 
characti'ristics  which  are  of  iihysical  interest  and  follow  from 
the  analysis  presented  herein, 

SYMBOLS 


coefficient  of  r"  in  velocity  profile 

(etp  (15)1 

ll. 

given  by  equation  (4.St 

lij 

constant  average  value  of  n. 

b 

positive  constant  used  in  reb-rence  2 

b„ 

coefficient  of  7"  in  stagnation-enthalpy 

profile  (e(|.  (Hi)) 

6. 

coefficient  in  thermal  profile  not  deter¬ 

mined  in  advance  by  boundary 
conditions 

A, 

constant  average  value  of  Ai 

r 

proportionality  factor  in  temperature- 

viscosity  relation  (('((s.  (11)  ami  (7)) 

r 

constant  average  value  of  (' 

T\-=  f‘f,/£ 

Jo 

average  skin-friction  coefficient  for 

length  L  (erp  (55)) 

<'ll 

local  skin-friction  coefficient  (e(|.  (.5li)) 

Cp.Cr 

specific  heats  at  constant  pressure  and 

constant  volnme.  respectively 

integrals  defined  by  e((nations  (i  1) 

/<■, 

constant  average  value  of  f’, 

Fu 

constant  replacing /’j  for  determination 

of  separation  iioint 

f  r ■  f Cl  (  ? ) .  f  r  4  (  ? )  ■  f  C.  ( 1 ) 

(larameters  defined  by  e(ptations  (51). 

(541)).  and  (4l)b) 

H 

stagnation  enthalpy,  (i/-/2)  +  Cp7’ 

h 

ratio  of  stagnation  ('iithalpy  at  wall  lo 

stagnation  enthaliiy  at  outer  etige  of 
boundary  layer,  //„///, (^);  for  Fr=  1. 
li=-TJTe  (cf.  also  eti.  (25)) 

be 

value  of  h  for  zero  heat  transfer  at  wall 

k 

L 

coefficient  of  heat  conductivity 
I'haracteristic  length 

M 

Mach  number 

III  .j.l 

constants  defined  liy  equations  (551)) 

\ii 

Xusselt  number 

Fr 

f^randtl  number,  licjk 

<] 

local  rate  of  heat  transfer  at  wall 

Hl 

Keynolds  number  based  on  L,p„u.^leln„ 

s 

r 

T, 

t 

II.  r 

j-.i! 

y 

V 

M 

£ 

s 

P 

T 

“rl 

yx 

Subscripts: 

a 

h 


<) 


1 


minimum  criliiid  J{c\ holds  iiuinhij’ 
based  on  conditions  at  point  h  nn- 
niciliatcly  behind  shock  uavc  at 
Iciulin'T  I'doc  of  airfoil,  p.. 

minimum  critical  Ucynolds  nninlici 
laised  on  remo!"  frcc->l ream  comli- 
lions  in  su|)crsonic  How  o\i'i  thin 
airfoil 

ratio  of  lociil  skin  friction  to  .\us-cli 
number  defined  in  c(|uation  itisi 
Snthci-huid  constant:  .S'  'Jlfr  K  for 
!iir  (cf.  cq.  (7) ; 
absolute  temperature 
equilibrium  wall  temperature  for  zero 
heat  transfer 

I  ransforma!  ion  \'a)i;i))lc.  ilclincd  lo 
eipiation  (Si 

vclocitx  components  in  r-  iind  //-dnec- 
tions.  respectively 

coordinates  parallel  and  normal  lo 
surface.  resf)ccti\ cly 
constiint  defined  by  equation  (:!4bi 
rat io of  specific  heal s. (•;,  r, : 7  1 .  f  for  air 

boundary-layer  thicknesses  in  .) 7  and 
.)'/  phines.  respectively 
recovery  factor  (c(|.  l4  I  )  1 

solution  for  lo  be  used  in  dctcr- 

minine- separation  point  (c(|.  iiHii 
coellieicnl  of  viscosity 
dimensionless  distance  alone  wall,  r  I, 
mass  density 

dimensionless  \iiriablc,  /  5, 
constant  dclined  by  eijiiation  l.■):^l 
constant  rephicine  v'l  i"  determinine 
sc|)aration  |)oint 

region  at  which  ad  verse  pressure  e|;i,li- 
ent  starts 

vjilue  at  point  outside  of  boundary 
layer  inmiedialely  behind  shock  wave 
at  Icadinj;  edee  of  su|)ersonic  airfoil 
value  at  wall 

value  used  for  determining  separation 
point 

vahu'  at  si-parat ion  |)oint 
local  value  at  outer  edec  of  boundary 
layer 

value  at  suitable  reference  |)oinl  out¬ 
side  boundary  layi'r;  in  numerical 
('xam|)les,  denotes  value  in  undis- 
tnrbed  (remote)  free  strt'am 


A  ])rime  denotes  difl’erenliation  with  respect  lo 

COMPARISON  OP  INTEGRAL  METHODS  FOR 
LAMINAR-BOCNDARY-LAYER  ANALYSIS 


Since  the  development  in  1921  of  the  boundary-layer 
momentum  integral  ecpiation  by  \'on  Karimin  (ref.  .'15)  tmd 
its  first  ai)pliealion  by  Pohlhansen  (ref.  Illi),  the  Karman- 
Poblhattsen  method  has  probably  been  the  most  widely 
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applied  and  fruitful  of  the  approxiinale  luetliods  used  for 
theoretieal  aiiahses  of  boundary  layers. 

The  Karnian  intefjral  e([uation  ean  be  lepirdt'd  physically 
as  a  moinentuin  balance  over  a  fluid  element  extendin*;  across 
the  entire  boundary-layer  thickness.  Matheniatically,  the 
equation  can  be  regarde<l  as  an  integration  of  tbe  original 
monientuni  partial  differential  eciuation  ov(‘r  the  boundary- 
layer  thickness.  The  advantage  of  this  intc'gral  e<|uation 
for  theoretical  calcidations  is  that  if  certain  <lefinite  forms 
are  assumed  for  the  velocity  profiles  a-  functions  of  the  nor¬ 
mal  distance  from  the  surface  then  an  ordinary  differential 
equation  is  obtained  with  axial  distance  along  the  surface 
as  indi'pendent  variable  and  essentially  the  boundary-layer 
thickness  as  the  unknown. 

The  Karman-Pohlhausen  method  in  its  original  form  is 
based  on  tin'  ust'  of  fourth-degree  velocity  pi'ofilcs  satisfying 
certain  conditions  at  the  wall  and  at  the  outer  edge  of  the 
boundary  layer.  By  means  of  this  |)articular  method  a 
considerable  variety  of  useful  results  for  laminar  boumlary 
layers  has  been  obtained,  even  for  cases  (d'  a  normal  mass 
flow  (fluid  suction  or  injection)  at  the  wall  witli  or  without 
lu'at  transfer  atid  pressure  gradii'iit  (cf..  e.g.,  refs.  .'{().  dT  to 
d9,  and  in  to  22).  It  has  been  found,  however,  that  this 
method  has  at  least  two  distinct  disadvantages  in  practical 
cases.  It  fails  to  predict  accurately  the  se|)aration  point 
in  an  adverse  pressure  gradient,  and  it  often  does  not  yield 
sufficiently  accurati'  results  for  derivatives  of  the  |)rofiles 
for  use  in  laminar-boundarv-layer-stabilitv  calculations.  In 
view  of  such  limitations,  various  refinements  in  the  Ixarman- 
Pohlhausen  method  have  been  made,  and  a  nutnber  of  what 
appeared  to  be  the  most  important  tyjtes  of  refinements 
were  studied  and  compare<l  in  reference  1, 

REFINEMENTS  OF  KARMtN-POHLHACSEN  METHOD 

In  discussing  refinements  of  the  Karman-Pohlhausen 
method,  it  should  be  first  observed  that  the  Karman  mo- 
nuMitum  integral  equation  is  not  actually  e(|uivalcnt  to  the 
original  partial  differential  eriuation.  It  is,  in  fact,  essen¬ 
tially  oidy  an  average  of  this  equation  over  the  boundary- 
layer  thickness.  Tims,  any  solution  of  the  partial  differ¬ 
ential  eciuation  will  necessarily  satisfy  the  momentum  inti'gral 
ecpiation  but  not  vice-versa.  This  iiasic  limitation  of  tbe 
integral  equation  is.  however,  to  some  extent  overcome  in 
tlie  Karman-Pohlhausen  method  by  the  fact  that  the  velocity 
profiles  which  are  assumed  in  this  ecpiation  are  not  chosen 
quite  arbitrarily  but  are  chosen  as  well-behav(>d  functions 
(namely,  fourth-degree  polynomials)  satisfying  the  boumlary 
conditions  and  certain  additional  conditions  which  an  exact 
solution  of  the  governing  partial  differential  ('((uations 
wotdd  necessarily  satisfy. 

There  are  two  main  types  of  methods  of  refining  the 
Karman-Pohlhausen  method.  One  method  consists  in  using 
integral  cfiuations  in  addition  to  the  Karman  momentum 
integral  obtained  by  multiplying  the  original  momentum 
partial  differential  equation  by  the  axial  velocity  v  or  powei-s 
of  II  (e.  g.,  refs.  40,  41,  and  29)  or  by  the  normal  distance  i/ 
or  |)owers  of  y  (ref.  42)  and  then  by  inti'grating  the  resulting 
equations  over  the  boundary-layer  thickness.  In  this  type 
of  method,  additional  unknown  parameters  as  finictinns  of 


tlie  axial  distance  j-  are  iiil |■o^bl<■ed  into  the  as--unic<l  \  e'ocil  \ 
profih-s,  and  tlicsi'  are  determined  l)y  the  aihlitional  resulting 
ordinary  differential  {‘(juations.  In  most  actual  applii  ations. 
oidy  one  integial  e<(uation  in  addition  to  the  Kiinmin 
monnmtum  integral  equation  is  introdueed,  and,  hence,  onl\ 
two  ordinary  differential  e(|uations  for  two  parameters 
result,  .'such  methods,  in  fact,  are  therefore  sometimes 
i-alled  “two-parameter  "  methods  .\  detaileil  discussion  of 
such  methods  is  given  in  reference  I. 

The  .second  main  type  of  refinement  of  the  Karman- 
Pohlhausen  method  is  the  use  of  onl\  the  Karman  integral 
(•((nation,  but  in  conjunction  with  profiles  cd'  higher  dcLriee 
than  the  fourth,  satisfying  additional  conditions  at  the  wall 
and  at  the  boundarv-laver  edge  which  an  exact  solution  of 
the  (tartial  differential  eipiations  woul  I  necessarily  satisfy 
I n  tnost  applications  of  t his  type,  velocity  piadiles  of  t  he  -ixt  h 
degree  (refs.  42.  44.  1  to  .2.  (i.  7.  2.2.  and  24i  are  used.  How¬ 
ever,  velocity  profiles  of  higher  degree  than  the  sixth  lia\e 
also  been  used  (refs,  41  and  22).  .Seven I h -degree  \cloeil\ 
profiles  have  been  found  pait icularly  suitable  for  calculation 
of  the  separation  point  in  an  adverse  pressui'e  gradient 
(refs.  45.  2,  and  7),  ( )ne-|)arameter  methods  with  \elocil\ 

profih's  of  higher  than  fourth  degree  are  discussed  in  some 
detail  in  reference  1 . 

COMPARISON  OF  METHODS 

In  view  of  the  variety  of  specific  means  of  refining  thi' 
Kiirman-Pohlhatiscn  method,  a  theoretical  investigation  of 
the  relative  merits  of  these  methods  was  tnadc  in  reference  1 . 
The  methods  weie  compared  on  the  basis  of  both  accuracy 
and  ('as(>  of  computation.  d'ln'  method  id'  comparison  wds  a 
posteriori.  A  relatively  sitnple  flow,  namely,  the  incom¬ 
pressible  and  com|)ressiblc  flow  for  a  Prandtl  number  of  unity 
iti  a  zcio  pressure  gradient  over  a  surface  at  a  unifortn  tem- 
I)erature.  was  calculated  on  the  basis  of  a  number  of  the 
foregoing  methods,  and  the  results  were  compared  with  the 
accurate  method  of  analysis  of  refeix  iice  24  for  flow  over  a 
flat  |)httc.  The  two-parameter  methods  considered  were 
based  on  (iti  addition  to  the  K.trmati  momentum  integiali 
the  integral  of  the  momentum  partial  differential  ('({nation 
multiplied  by  n  in  conjunction  with  fourth-  and  fifth- 
(l(‘grec  velocity  (irofiles.  The  one-paratneter  methods  were, 
of  course,  based  on  the  Karman  momentum  integral  e(pia- 
tion  and  were  a])i)li('(l  in  conjunction  with  fourth-  (Karman- 
Pohlbatisen  method),  fifth-,  atid  sixth-degree  velocity  |)ro- 
files.  Th('  comparisotis  wen-  tnadc,  in  ()artictihtr,  on  the 
basis  of  calcitlated  skin-friction  and  beat-transfer  coefficients, 
first  and  sccotid  derivatives  of  the  profiles  throughout  the 
boundary-layer  thickness,  and  minimum  critical  Ki'vtiolds 
numbers  for  laminar-boundary-laver  instability  accorditig 
to  the  criteria  of  Lin  and  la'cs  (refs,  4li  and  .Hi.  It  is  well 
known  that  the  latter  criteria  arc  sensitive  to  first  and  s('con(l 
derivatives  of  the  |)rofilcs. 

The  botitidary  conditiotis  satisfied  by  tbe  various  (uofiles 
as  well  as  the  detailed  conqtarison  of  the  restilts  of  the  various 
methods  can  be  fotitid  in  referetice  1.  Tbe  results,  in  brief, 
indicated  that  skiti-friction  and  heat-transfer  coeflicients 
were  pn'dicted  with  stibstatitially  satisfactory  acctiracy  by  all 
of  the  methods.  Moreover,  the  overall  profile  shaitcs  ob- 
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taiiied  by  all  of  the  methods  were  <iiialitatively  eorreet.  How¬ 
ever,  quatititative  tlifferenees  in  the  first  and  espeeially  second 
derivatives  of  the  profiles  were  obtained,  with  correspondinf; 
differences  in  the  calculated  values  of  the  minimum  <-ritical 
Reynolds  numbers.  It  was  concluded  that,  on  the  whole,  the 
one-parameter  method  with  si.xth-degree  profiles  gave  the 
most  accurate  results  for  the  profile  derivatives,  as  well  as  for 
the  minimum  critical  Reynolds  numbers.  In  reference  1  the 
stability  calculations  were  carried  out  for  the  case  of  zero 
heat  transfer  at  the  wall.  Subsequent  calculations  (ref.  3) 
indicated  that  reliable  results  for  stability  caleulations  by  the 
one-parameter  si.xth-dcgree-profile  method  are  obtainable 
also  for  the  case  of  heat  transfer  at  th(>  surface  of  the  flat 
plate. 

In  addition  to  being  capable  of  yielding  results  of  adeipiate 
accuracy,  it  is  usually  also  quite  desirable  that  a  method  of 
calculation  be  simple.  In  this  connection,  it  must  be  ob¬ 
served  that  the  one-parameter  methods,  in  general,  involve 
considerably  simpler  calculations  than  the  two-parameter 
methods.  This  ailvantage  of  the  one-parameter  methods 
may  not  b('  very  pronounced  in  the  case  of  flow  over  a  flat 
plate:  however,  it  becomes  quite  pronounced  for  th(>  general 
case  of  flow  in  a  pressure  grailient  with  heat  transfer.  In  this 
case,  the  thermal-energy  partial  diffen'iitial  eipiation  must  lx* 
integrated  to  yield  an  integral  equation  in  addition  to  the 
momentum  integral  equation.  Consequently,  there  will  be 
at  least  two  parameters  to  determine.  If,  however,  both  the 
momentum  and  tin-  thermal-energy  partial  differential 
('(luations  are  multiplied,  for  example,  by  ii  and  integrated 
over  the  boundary-layer  thickness,  tlp'ii  a  total  of  four 
ordinary  difl'erential  ('quations  in  four  unknown  parametei-s 
will  be  obtained.  Thus,  the  so-called  two-parameter  method 
wouWl  in  this  case  really  become  a  four-parameter  method. 
iThe  one-parameter  method  in  this  general  case  similarly 
becomes  a  two-parameter  method.)  It  is  noteworthy,  iti 
fact,  that  in  any  of  the  foregoing  applications  of  the  two- 
parameter  method  oidy  the  less  general  cases  of  zero  pressure 
gradient,  or  pressure  gradient  with  zero  heat  transfer  at  the 
wall,  have  been  treated.  If  it  is  desired  to  develoj)  a  unifieit 
method  to  be  applicable  in  the  more  general  as  well  as  in  the 
simpler  cases,  then  this  would  have  to  be  considered  a  dis- 
atlvantage  of  the  two-parameter  methods. 

In  view  of  the  foregoing  results  and  considerations,  it  was 
concluded  in  reference  1  that  the  most  [iromising  integral 
method  for  laminar-boundary-layer  study  appi-ared  to  be 
that  based  on  the  Karmi'm  integral  equation,  in  conjunction 
with  sixtlwlegree  velocity  profiles.  This  is  essentially  the 
method  of  analysis  to  be  applied  in  the  present  report.  It 
should  be  observed,  however,  that  cases  exist  in  which  profiles 
of  other  degrees  are  preferable.  Iti  particular,  stagnation 
flows  are  more  satisfactorily  treated  by  fourth-ilegree  profiles 
(ref.  2).  while  the  separation  point  in  an  adverae  pressure 
gradient  appears  to  be  determined  more  accurately  by 
scventh-ilegree  profiles  (ref.  2).  The  latter  case  will  be 
treated  in  some  fletail  in  the  present  report. 

It  may  be  asked  why,  in  the  one-parameter  method, 
profiles  of  higher  degree  than  the  sixth  were  not  considered 
in  the  comparison  study  of  reference  1.  The  reason  is  that 
the  sixth-<legree  profiles,  as  distinguished  from  fourth- 


degree  profiles,  an'  cho.scn  to  .satisfy  an  additional  condition 
at  the  wall  (as  well  as  the  outi'r  boundary-lax cr  edge).  'I'liis 
condition  is  obtained  by  ilifferentiatiug  the  jiartial  difb'rci.tial 
momentum  equation  with  i-espect  to  r.  If  a  velocity  profile 
of  higher  degree  than  the  sixth  is  assumed,  then  the  onix 
means  of  obtaining  a  further  eoialition  at  the  xvall  xxliieh 
woidd  be  sati.sfic<l  by  an  exact  solution  of  tlie  partial  differ¬ 
ential  equations  is  to  differentiate  the  momentum  partial 
xlifFcrential  equation  txvice  xvith  res|)eel  to  t  (or  //  foi-  in¬ 
compressible  floxvs).  and  then  take  values  at  the  xvall.  This, 
however,  xxill  be  found  to  yiehl  a  condition  involving  |)ai  tial 
derivatives  xvith  n'spect  to  r  such  as  [d-iiij.i/)  dj-  and 

this  I'ondition  then  becomes  essenti’dly  an  addil  ion.il 
ordinary  differential  equation.  Since  the  sixth-degree  pro¬ 
files  hax'e  apparently  led  to  satisfactory  residts.  it  has  not 
seemcxl  worthxvhile  to  introduce  such  eomplieat ions  into  the 
analysis  by  iisitig  bigher  dcgri'c  velocity  profiles.’  It  is 
notexvorthy,  in  this  regard,  that  'dthotigh  |)olx  tiomials  of  as 
high  a  degrt'e  as  the  elex'cnth  xvefe  apiilied  in  referenee  41. 
they  .satisfied  only  the  same  eotiditiotis  at  the  xxall  as  the 
sixth-degree  jirofiles  to  be  ti.scd  in  the  itresent  re[)ort . 

BA.SIU  EtH'ATIOXS 

The  folloxving  e(piations  describe  the  steady,  txxo- 
dimcnsional,  laminar-botmdary -layer  lloxv  of  a  cDinpressihh' 
gas  along  a  slightly  etirved  xvall: 


E<iuations  (I),  (2),  and  (4)  are  the  momentum,  continuity, 
and  energy  e(|uations,  respeetively.  Erpiatioti  (3)  folloxvs 
from  the  ideal-gas  laxv  and  the  assumption  that  the  |)resstire 
is  constant  across  tin*  botitidarv-lavi'r  thiektiess.  It  xxill  he 
a.-ssumed  here  that  the  spi'cifie  heals  Cp  and  Cc  as  xvell  as  the 
Prandtl  number  I‘r  are  conslatits.  By  multi|)lyitig  eqtiatioti 
(I)  by  II  and  adding  the  restilting  x'tpiation  to  ecptalion  (4). 
the  folloxving  form  of  the  etiergy  erittation  is  obtained  for  a 
constant  Prandtl  number: 


lr[pu 


~l‘r)  (a-: 


It  will  be  assumed  in  the  present  utialysis  that  the  viseosity- 
temperature  relation  can  be  a|)proximated  in  the  form  (ef. 


refs.  34  and  15) 


(«) 


‘  Whvn  sovpmh-iJoBHf  vulnrUy  pntfllps  uru  list*'!  litTf  f«*r  calculation  of  the  suparaiioti  I'oini 
in  an  atlvcryc  pressun*  .gratliont,  the  a<l*liflonal  cttinlltion  satisflutl  at  the  waM  is  chosen  to  be 
exactly  vaH*!  only  at  the  separation  p*)lnt. 
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(T) 

Ii  is  coiivcniiMil,  in  this  (•oiii|)n'ssil)lc-fl>iw  aiiidysis,  lo 
r(>pla<‘i'  iho  iionuiil  distaiicc  cDonliiiiitc  1/  i)y  llu'  Dcinxliiil/yii 
Vid’iahh'  1  (l(‘fiiii'(l  l)v 


By  iiitcfrratiiii;  o(HiatioMs  (1)  and  (n)  witli  respect  to  i  over 
the  houndary-layer  thickness  /  -(I  to  t—&,  and  iisin<;  the 
honndary  conditions  1/ -  r  ()  at  t  ()  tofretlier  with  smooth 
traitsition  of  the  velocity  and  t(‘in])('ratitri'  jtroliles  to  their 
local  main-stream  valni's,  the  followini;  integrodilferenlial 
e(|uations  are  obtained: 

{t\  'll +"|-  [^■.+(>+^._,  ' 

:('(p,  P,)(»,  lh)n\  "dj  ('.*) 

-t'tp,  p,)(»:,  )/,)(7’|  7,  )U  //|)J  (1(1) 


t\^  (u  )iy)  [1  —(u  a,)]  dr 

fV— J  \(II  Ih)  —  {n  'ti)'’|'/r 

('/  '/i)  (I -(//7/i)|'/t 

x=(a,,X)-(p., 


an«l  21  ).*  Moreover,  as  explained  in  referenee  (i.  llie  use  of  a 
sini;le  houndary-layer  t  hiekness  does  not  lua  essarily  nuposi' 
any  undite  ri-st riet ions  011  the  thermal  prolih's.  sinee  the 
'  hitter  have  here  heeii  permitted  to  contain  an  ailditioiial  eo- 
efheient  not  determineil  in  advanei’  l)\  the  lannidarx  condi¬ 
tions,  This  eoeflieient.  to  he  taken  here  as  /i|.  replaces  the 
thermal  honiidarx -layer  thickness  as  the  second  unknovx  ii  to 
he  determined  hy  equat  ions  ((t)  and  (  HI )-  A  sinoh'  hoiindai  >  - 
layer  thickness  has  also  heen  used  in  i(‘fer<-ni  e  22. 

!  Kquations  {!!)  atid  (lOl  can  he  converted  into  ordinary 
dili'erential  equations  hy  assuming  the  velocity  and  sta;_nia- 
I  tion  enthalpy  as  definite  functions  ol  the  normal  disiaiiee 
lariiihJe  r.  For  this  /lurpose.  as  ev|)lained  in  tlie  section 
I  "( 'om|)arison  of  Intcirral  Methods  for  l,aminar-Boundar\ - 
j  Layer  Analysis,"  the  \eloeity  jiroliles  will  he  chosen  a- 
si.xtli-deirree  polynomials.  'I'lie  sta<2nat  ioii-eiit  halpy  pndiles 
will  similarly  he  chosen  as  polynomials  hut  of  one  de;_uee 
hi<rher.  namely,  seventh  decree. 

'riius,  it  will  he  assumed  that 

a  H I  (/,.  r"  '  I  ■< 


II  IF 

n  --11 

'I’he  following  honndary  l•onditions  iniisi  he  satisfied: 
At  r^  I). 

a  r  t)  'I 


Here  X(|)  is  essentially  the  nonditiu'iisional  s(|uared  houtidarv- 
layer-thiektiess  jiarameter  iti  tl.e  st  iilane. 

The  quatitities  p,'p,  and  .17,  in  ecptations  (tl)  and  (HI)  are 
related  to  which  is  a  funetioii  of  ^  preserihed  hy  the 

potential  flow  ahout  the  hody  in  (piestion.  Thus,  in  aei'ord- 
anee  with  tin'  usual  isentro|)ie-(low  relations, 

I  I 

PiP.  ==(7’,  '7’J>->  =={1^(7- 1)  (.1/7^2)11 -(»,/»  (L!) 

.U,=  (a,/aJ.\/„(7V7'J-''-  (14) 

In  deriving  equations  (!))  atid  (10)  ti  single  honndary  thiek- 
iK'ss  has  heen  assumed.  This  is  an  alternative  to  the  inlro- 
ditetioti  of  two  houndarv-layer  thicknesses,  namely,  a  ve¬ 
locity,  or  dynamical,  and  a  stagnation-enthalpy,  or  thermal, 
houndary-layer  thiektiess  (ef.  refs,  10  to  21  and  2d).  The 
assumjition  of  a  sitigle  houndarv-layer  thiektu'ss  appears 
feasible  for  fluids  with  I’randtl  tiumhers  close  to  unity,  since 
iti  that  ease  atialyses  involving  both  a  ilynamieal  and  a 
thermal  houndarv-layer  thiektiess  tisually  imjily  that  both 
thicknesses  are  approximately  etpial  (see,  e.  g.,  refs.  dO.  20, 


II  IF  f>(i) 


where  /(({i  is  eonsith'ied  as  a  presinihed  function. 

At  7-  I , 

,j,i,-IIIF  ^  1 

X  I  IS) 

dll  0  -—dll  dr  0  J 

Iti  addition  to  these  conditions,  the  following  . . litions  will 

also  hi'  satisfied  (ef,  ref.  (i): 

At  7^=0. 

(VI'-,  TJ  ^  ..{wild  —  X(p,  p  ,  )/i  |1  :  (7-l).U,-  2|((/,  u  ,  )' 

07“ 


d'  ,  ,  dHu  II,)  dill  IF) 


dF  07 


i  (7-f).Ur'2j  ",(////,)  -  (I -/V)(7-i).\/r 

0  7" 


ro(  II  '/i)T 
‘  L  dr  _ 


{ii,'ii  j\l‘r{p,  pj 


d(l-/V)(7-l).\/ 

At  7=-  1, 


7'r’eerr,7,.{^";.(////,.- 

,  ,2,,, 

Or  07“  J 


(iiiii,)  --^^.^  ("/Mi)--^^2  (^77/1)-^^^,,  (IFIF)^d  (2d) 

*  In  (»f  titiw  ni'ttr  i\  stiiynniioii  iHiini.  Iiowpvrr,  it  is  iim  lostini;  in  iiulf  lh;»i  i  vri)  lor 

SI  I’nuKltl  tmnilxT  nf  I  tho  ihfrnv.il  lonintijiry-liiyor  thick tu'ss  in  the  n  phinc  nm;.  h.  :ir<nin<l 
2/)  iMTtvnt  fircalcr  than  ih'-  -lytminical  htiuntlury-laycr  thickness  (ref  •.‘tn. 
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( 'olitlil  itiiis  il!ti  til  ll’2i  I'lillinv  iVoin  l■l)Ullli(llls  II  i  to  i4i  iiiul 
iliH'i'i'i'til iatiiiii  Ilf  I'licli  (if  (■(|iialiiiii>  ili  anil  ifi  (incc  with 
l•l■s|l(‘('t  til  r.  takin>:  values  at  tlie  wall  ami  takiii"  enmlitidiis 
(  1 7 1  in  til  aeeiinnt .  ( 'iinilit  inns  i  J.'!  >  fnllnw  frnin  ililferent  iai  inn 

nf  e(|natinns  (1)  ami  (4i  nnee  ami  twiee  witli  res|)eet  tn  t. 
takinn  values  at  tlu'  Ineal  miter  linnmiatv-la yer  eiljre  ami 
takine  enmlitinn^  (ISi  intnaeennnt. 

After  till'  enellii  ielits  ii„  ami  /  have  lieeii  ilelertnineil  frnm 
enmlitinn>  (  17i  tn  r.‘:i  i  in  tertns  nf  o,  inr  Xi  ami  h,.  e(|iiatinns 
ill  I  ami  (It)'  lieenme  nnlinarv  ililferent  ial  e(|natinns  tn 
ileterinine  X'Ji  anil  Fur  any  iiiveii  ease,  the  llnw 

nntsiile  the  hnnmiary  layei',  as  ilelineil  liy  n,  n  t^i  ami 
M  .  is  ennsiilereil  as  |ireseriheil.  .\lnrenver.  the  tempera- 
tnre  ilist rilint inn  alnnii:  the  snrfaee.  as  iletineil  liy  /I'ti.  is  alsn 
ennsiilereil  as  preseriheil  here. 

'The  temperature  pinliles  are  relateil.  in  ireneral.  tn  the 
stajiliatinn  enthalpy  ami  the  velneily  pinliles  in  aeenrilanee 
w  ilh  the  relatiiin 


Ki|Uatinn  !'J4i  Inllnws  frnm  the  ilelinitinn  nf  the  >tairnatinn 
enthalpy  //.  Frnm  etptatinn  il'4'  it  fiillnws  that  the  wall 
tetnperatnre  (list I'ihiit inn  is  relateil  tn  the  sti  j:natinn- 
eiithalpy  ilistrihntinn  /m^i  at  the  wall,  in  aeenrilanee  with 
t  he  relal  inn 

hii){\+\,'  Mr) 


I'rnliles  in  the  .li  plane  can.  if  desirei.  he  I  ransfnrmeil  iiitn 
pinliles  in  the  |)hysieal  .m/  plane  by  ileterminin^  //  as  a 
funetinn  of  r  ami  J  in  aeenrilanee  with  eipiatinn  iSi'. 
Kipiatinn  (Ni  ean  he  expressed,  in  General,  in  the  fnllnwiiur 
nondimensinnal  form  : 


0 


p  r  T.rlr 
I  '  iTTr 'It 


tL’t'iai 


where  i  is  t  he  physieal  hnumlary-layer  t  hiek  ness  determined 
by 

yx|  '  l-’libi 


attaek,  i-  dev  elnpi  d  'rin-  I’l  aiidll  iinmbri  i~  c  iii I'li li  i  i  i  1 
arbitraiA  but  ennstalil  and  nf  ihe  mdei  nf  m.’e.:n 1 1  ud>'  nl 
unitv.  while  ihe  uall  lemperalnie  n!a\  \ai  i  alnny  ihe  linu  ; 
1  hat  i>.  i .  y  :  nr  A  A  A 

(iKSKKYI.  sni.l  nns 

Fnr  l|n«  wilhniil  a  pre — me  yiadieiil,  ihal  i'.  n  i' 
while  ,l/i  .1/  .  /|  /  .and  >n  fnrih.lhr  in  l  eul  nd  ill .  i  e,  u  i.i  I 

eipiat  inns  i  ‘I '  and  MU’  beennie 


The  si\|  h-dee|ee  \e|neily  pin.ile  ~;i  I  l-l  \  lie.;  bniniil.iri 
enm lit  inns  (17’  t  n  ’  I’n  ni  I  h  i-  l■.■|~r  i- 


i; 

II I 


'it  Ur 


The  M-velit  h-denree  -I  aei  la  I  ini  i-ei  1 1  ha  I  p\  pinlile  -.■i  t  mI  \  1 1 1_' 
these  enndillniis  in  ihl-  ea-e  I'  L:l\eii  b\  iiiua'lnn  lU 
w  here  '■ 


A„  //„//  A-t- 

A,  •><;_ 
h;  -.’fl'  X 

/),  :;.'iM  A I  2UM’  sM' \  2uA. 

A,  Nil  I  A’  lut;.  IL'M'  X  l.'.A 

A„  7U(  1  III  ::uf«',  Nfi'.x  :;uA, 

!>:  2U1  I  A  -  St;  ■  ■>(;  X  ■  IIIA, 


>  UU) 


J 


w  here 


I  I  /'(  I  I'l  I  1.1/ 


:il 


In  the  sueeeedinir  twn  seetintis.  it  will  be  shnwn  how  the 
eipiatinns  thus  far  developed  ean  he  used  In  lead  In  a  simple 
and  usually  sudieienlly  aeeiirale  methnd  nf  ealeiilatin^  the 
laminar-bnuiidarv-layer  eharaeteristies  fnr  eerlaiii  "eneijil 
ly|)es  nf  flows, 

FLOW  WITHOIT  AN  AXIAL  PKKSS14JK  (tKAOIENT  WITH 
AUBITRARY  (ONSTANTi  PRANDTL  NTMIIER  AND 
VARIABLE  WALL  TE.MPERATIRE 


and  />_■  is  thus  a  preseriheil  eou'lanl.  w  hile  (I.  is.  ni  general,  a 
<riven  funetinn  nf  i;.  Frnm  ei|ualinn>  Ml  m  eniijmielinn 
with  e(pialinns  (2'.>i  and  Ultli  it  is  Inimd  that 


/•| 


•IS.') '.I.uil'.l  I 


Ai- 


!i.'):i 

LSI).  LSI) 


:102  s-Ji 

ii.lMMi  ''“'12.012 


a 

j 


I 


1 


In  this  seelinn,  based  on  the  e(iualinns  thus  far  derivi'd. 
a  simple  and  aeeurale  method  fnr  ealeiilatiii}?  the  lamiiiar- 
hnumlary-layer  eharaeteristies  of  the  (low  in  a  zero  pressure  | 
nradieiil.  such  as  the  llnw  over  a  flat  plate  at  zero  anetle  of  j 

J  'l  lii.s  \s  ti(t(  tlfci'.sarx  il  it  U  tlusiri  -1  t'l  'iptiTinirif  only  i>ropc!  tics  such  a-  Nkiii  frir'inri.  hr  it  [ 
(r  insf.r,  'iri'l  srpinition  point  which  licpcM'l  only  on  vnl'ics  tu  'ft.  wall.  ' 


With  /*’,  as  nj  Veil  bv  eipiat  inn  (22  ) .  I  he  sol nt  inn  nl  e(|ual  imi 
(27l  for  X(p  with  the  eouditinn  X(0)  -  I)  is 

X  4/',  /•',  Uhfi 

•  'I'h.'  syliihnf  ol  rofiTcnccs  li  tUbi  7  i  roplm'.-'l  tu  ri'  liy  !  In-  si  iHiimy  h:ii  tn>u  f  ipi  u  •  d”  mi. 
s>  nihi.l  h  inscl  also  in  refs,  'jii  ami 


s 


i;  I-,  pule  l.'l.'i 


\  A  nu\ \i.  \i>\  l>ui;\  i  iiMMi  r  i  l  i  I  ui;  Ai  l:u\ Al  I  II  " 


wlicli' 


Mol'i'i>\  IT,  willi  h\.  ti.  iiiiil  X  ii>  ^ivi'ii  liy  rij  Hill  iuiis  u{2i  iiiiil 
llic  solution  of  l■<llilllio^  (2S!,  wiili  the  eoiuiiiion  ol 
liiiili'  />!  ill  t  <1,  is 

X,  (\  ■’ I  i/J  (u-fai 

where 

a  I  1  2l  •  ili.s.')  2.4li:{)(  I  /'/  I 

fi'.it)  il'.Mil'J  s-Jliicil 'J.v.'Ml  /,)(rr,i  A'  cu  i-.'iii  ■ 

^  i;uh) 

iii.Aii/'/'  J!i,.‘i.')iii|':i  •_>!  \<\i" 

I !i.'>:i/'/' 2!i. . ■).■>() II r,  riA"  Moi  ii.iiii'.iu^ '  r,if;,)^ 

From  e(|iuilions  CUiii  iiml  (o4h|  with  !iu\  preseriheil 
I  einpi'l'ill  lire  ills!  rihillioil  Alii  llt  the  Willi,  A,  A|I{|  eiiii  he 
reii(lil.\  ileteriiiineil  hy  ii  simrle  i|il!i(li  iil  lire  (for  whieh.  e. 
iiiimi'iieiil  inteoriil ion  mny  he  nseil).  From  e<|niilioiis  l.■!(ll 
the  remililline  A„  eoellieieiils  cilli  liiell  he  foi’inl  its  fniieliolis 
of  'iml  the  veloeily  nml  si  iionii  I  ioii-eii  Jiiilpy  profiles  in 
the  ir  plillie  lire  llieli  ilelermilieil  i>y  eipiiilions  (2U)  iiml  Iflii. 
The  tempeninire  profiles  follow  from  eipmlioii  (24).  'I'he 
profilt'S  eiin  he  I  riiiisformeil  inlo  the  physienl  .ri/  phine  h> 
menus  of  1  ipiiilion  |M  or  (2tii. 

The  iivei'iioe  shin-friel ion  eoeflieieiil  for  the  leinrth  /.. 
neeordino  to  enuiuiuus  (2'.li  ami  will  he 


.'Slllisl  il  III  me  now  A  A, I  7',  ;ef.  ei|  J.'i  :  '  llllo  e(|n:lli.i!. 

i.’i4h)iin(l  also  -.iilisl  1 1 'll  me  expi'e'sion  ii;'.l  loi  A  in  leini'ul 
fu.  il  is  loimil .  a  il  h  ^  rephieei I  h\  ;in  ii \  i  riiiie  eo n-l  in 1 1  \  ;i I  le 
i\  for  simplieii  \ ,  thill  eipiii  I  iuii'  i:!4iii.  li.'l'.  niul  T  \ie|,! 
the  folloa  ine  I'Xpression  for  llie  \ns>ell  iinmher. 

-Vi/  il.2‘i7\ (  I  -  y.' )  i  ’  ‘  \  A'/  I  'h  H'ii 

u  here 

<!:.  '  [1  1 7',,  7’.  I  !2  '(iii,i):;n  i:,2.i;7:u/’(|  7'  /'.  t 

I4ii.02ii  l.')2,r>s7ii/'i  I  7',,  7',  i"i'  hh' 

For  liny  preseriheil  ilisi  rihni  ion  of  ihe  anil  i  einpei  ii  I  nii'  ii- 
•riven  hy  /’„  7,(^1.  the  loi  ill  .\nssell  mimher  i  iin  In  re;iclil\ 
ohliiineil  from  eipnilions  I4lliii  iiml  i4ilhi  Tin  iieiii.il 
heiil-lninsfer  rule  ill  I  he  anil  i  in  niiils  sin  h  iis  Ui  n  per  'eeuml 
per  sqn lire  fool  i  eiin  1  In  n  he  ohi  iiiiied  imiiii'iliiilel  \  h\  ~iiK  iii” 
fori/  Ml  lerms  of  .\’|/  m  iieeoi'ihinee  anil  eipnilioli  '27 

The  eipiilihrinm  inliiiha  I  ie  anil  lem|ieintnre  7’  iniuiilinu 

to  eipnilions  i;pn.  ,21!.  nml  '2.'>'  is 

/:  ■/'[ I " 

a  here  i),  klioa  li  iis  llie  lemperntnre  reeo\  i  i  \  fnelor.  i-  loinnl 
to  he 

r;  i  II  272 1  I  /7  I  12 

All  exiiel  iiniilssis  le.  e,.  ref.  ,20)  slio.'  s  thii!  fo|-  line  a  illioiil  n 
pressure  ei'mlieiil  over  iin  impermeiihle  >nrfiiee  n  i  i'IA  eomi 
iipproximntion  for  i;  Ilo  he  iletioted  here  iis  q,,  1  is 


where 


I  I  i,u  0'/ 

1.  1 1  2 )  p  II  ■ 


1.222 

nA’. 


A  <  \ 


i2.'» 


(2li) 


The  .N  us.sell  niiiiiher,  a  hieh  is  n  nondimeiisioiiiil  meiisnre 
of  the  rule  of  heiil  Iriinsfer  ill  the  anil,  eiiii  in  this  eiise  he 
defined  ns 

.V// =  I//, /■  (  /',  /’„)  (27 1 

ahere 

1/  ( /■  d /' O// h,  i2Sl 


i),..  X /'/  12 

For  /’/•  0.72,  ahieh  is  esseiitiiillx  the  vnliie  for  iiir.  eipnilioii 

(42)  >  ields  1)  0.024  insteiid  of  the  neeni'iile  \nlne  rf.  -  0,s4.‘i 

(ref.  24i.  This  iniieeiiriiey  in  the  vnlne  of  q  implied  h\  ilie 
eipiiitions  used  here,  hoaever,  will  not  neeessnrilv  iill'eel  the 
iieenrnev  of  eipiiifion  l40n)  for  the  Xnsselt  nnniher.  sinee  the 
ileriviition  of  this  eijniition  aiis  netnnlly  miide  imlepemhmt I \ 
of  the  piirlienliir  vnlne  of  ti,  !i.  e..  imlepemlently  of  eip  '21 
mill,  lienee,  of  q:  This  is  further  verified  hy  the  iiirreerneiil 
ohtnineil  aitli  eertnin  exnet  solutions,  to  he  diseiissed  -nh- 
seipieiitly.  Tims,  in  npplyin>r  eipintions  (40n)  and  (27  i  for 
the  enlenlntion  of  hent-l rnnsfer  rntes.  the  iietnnl  value  of  /' . 


iind  /',  is  the  eipiilihrinm  anil  lemperntnre  for  zero  heiit 
I  rnnsfer.  that  is,  for  i/  (I.  It  is  ii])])ropriiite  in  the  deiermmn- 
lion  of  heiil-triinsfer  rntes  In  reitlnee  the  tem|)enitiire  pnrniii- 
eter  A  hy  the  |)nrnmeter  7',,  7',,  whieh  is  the  physienll.v 
sifrnifienni  temjierntnre  |)nrnnieler  in  hinli-speeil  flows  with 
hent  Irnnsfer.  Ati  e.xpression  for  ihe  Xnsselt  iiumher  in 
terms  of  7',,  T,  enn  he  ohtnineil  h\  first  finiliii<r  the  value  of 
A  (to  he  denoted  as  A ,)  for  zero  hent  transfer  at  the  wall.  By 
s.ihstitiilinjx  A  A,  into  eipialion  (24h),  assiimin;^  A,  eonstaiil. 
and  .set tine  O  (fm-  /.ero  heat  transfer,  or  A,  0)  the  foliow- 
ine  value  of  A,  in  terms  of  (1.  is  ohiained: 


ns  ilelermilieil  either  hy  experiment  or  h\  eipnit  ions  (4  I  !  and 
(42),  should  he  used. 

It  shinilil  he  ohserved  that  the  use  of  the  npiations  de¬ 
veloped  here  is  not  restrieted  to  any  partienhir  t.'pe  of  tem¬ 
perature  distrihiition  A(^)  1110110  the  aiill.  Thus,  it  is  not 
iieeessury ,  in  a\)pl\  the  method  of  eulenlalion  desevihed 
here,  tinit  the  temperature  distrihiition  he  expressed  ns  a 
polynomial  in  J  (unlike  ref.  24)  or  iis  a  poaer  of  t,  dinlike  ref. 
47).  Ill  the  speeial  ease,  hoaever.  in  ahieh  7'„  is  expressed 
as  a  polynotiiial  in  the  ealeiilat ions  indieateil  hy  the  present 
j  method,  ineliuliii”:  the  transformation  from  the  j!  plane  <0 


If  iu*vrrfhfl«*^s  tlpsiixHi,  :i  prat  aihI  \ n  y  siiiipli'  itican'i  nf  iiimiifv  iii-:  th<‘  pn  'on!  .npii- 


j  in  iiriln  virfimlly  lo  -liniinAfi'  th«‘  implit  ii  iiiscn-p;ui<*y  Imiwi  oti  r?  aii'l  »j.»  (|nr  /’r  ji.nr 

!  unity)  vviinlil  mrri'ly  to  mtiliiply  thi- valnr  nf  f/;' -v>  ui  vn  by  niirui.in  (dl  b\  J.  Kijiioi  i.itfs 
■  flUa.  and  ituf!-  u«a>!'l  n-itiain  uiicliatiui-'i  !•>  -aicii  a  tini  •  (■•ion. 


ANAI.VSIS  OK  LA.MIN'AIi  I'OM I’KKSSl lU.K  HOI  NDAUV  l.A^KK  WH  II  IIKA  T  I  lfANSl  KK  AM)  I'lil  SsrUK  (.IIADII  NI 


(I 


tli(>  plivsiciil  j’ji  plain',  <’an  In  cairiotl  out  (lircctlv  witlioiit  any 
ipiadraturi's  liy  iisinir  tlu'  results  jriven  in  tlie  appeiulix  an«l  i 
(ifrures  4  and  a  of  referi'iiee  (i.  j 

<’OMPARISON  WITH  KXAC  T  SOl.l'TION.S  , 

As  a  ehi'ck  on  the  aeeurae,\  of  the  results  ohtained  here, 
eomparison  has  heen  made  witli  eerlain  known  exaet  solii-  j 
tions. 

For  tile  speeial  ease  of  a  uniform  wall  temperature  and  a  j 
1‘randtl  nuiuher  of  imit\,  it  is  well  known  that  the  eiieifrv  1 
partial  tlill'erential  ecpiation  (4)  or  I.")')  reduces  exactly  to  a  i 
sineh'  (piadratie  relation  hetweeii  the  temperature  and  the 
veloeit  .  throughout  the  houndary  layer,  'Phis  relation  lan 
he  expressed  in  the  form  I 

ii  A  •  (1  - /-)(  ''  )  (44)  I 

By  puttine  A  -C’onstant  anti  /'/  I  into  etpiations  CtManil 
i:f4h)  the  solution  for  A,  as  aiven  hy  etpiation  Cfda)  is  founil 
to  he 

Ai^L'd-Ai  (4.')| 

Suhstitution  into  etpiations  (lid)  for  values  of  A„  and  eompari-  | 
'Oil  of  the  resulting  sta5;nation-enthalp\  proliles  with  the  ! 
M'loeity  proliles  Cj;))  then  show  that  relation  (44)  is  exactly  I 
'iitislied.  'Phils,  the  ('(piations  used  here  reduce  to  the  t'xaet 
intesrral  of  the  t'Uersry  (initial  tlill’erential  etpiatitm  in  this 
special  ease.  .Vs  alrt'atly  indicated  in  the  section  "(’omiiari- 
son  of  Integral  Metlunls  for  l.aminar-Bouiitlary-l.ayer  .Viialy- 
'is,"  it  has  heeii  ftmiul  iref.  1).  inoreox'er.  that  the  skin-  i 
friction  anti  heat-transfer  eoellieieiits  olitainetl  hy  the  (ireseiit 
iiiethtMl  in  this  ease  aoree  aliimst  exactly  with  tht>se  olitaiui'tl  ; 
hy  I  he  exact  inethotl  of  refereliee  (14.  'Phe  (ireseiit  inelhtal  ! 
Inis  also  heen  ftiiind  Iti  yield  results  of  satisfactory  aeeiiraey 
for  stiihilitv  ealeiilations  in  this  ease  (refs.  1  anti  If). 

'Po  cheek  I  he  resiills  tif  the  (ireselil  iiielhotl  for  the  limre 
ireiieral  ease  tif  /'r^^l  and  variahle  w  all  teni()eral  lire,  ealeiihi- 
lions  were  carried  out  for  the  ease 

7’„  7’,  -  l.l’.’)-  (I.N.'li  ■  (i.;i;!i" 

Phis  is  the  ease  ealeiilaleil  in  reference  14  hy  the  exact  i 
inelhotl  of  analysis  iheie.  'Phe  local  .Niissell  nuiuher  for 
this  ease  was  ealeiilaled  hy  iiieans  of  etpiations  (4(la)  anti 
(4(lhi.  In  atltlilion,  tein()erature  anti  velocity  (irofiles  were 
ealeiilaleil  hy  means  of  the  (iresent  etpialitms.  'Phe  ajrree-  , 
meiil  hetwi'cn  the  results  thus  ohiiiinetl  anti  (hose  in  refer- 
eiiee  .■i4  was  found  tt)  he  t|uite  close  (see  ref.  (i  for  iletails  tif 
ihe  ealeiilatitins  anil  resulls). 

Il  is  inlerestintr  to  note  that  hy  selling  0  anil  solviiifr 
ihe  result iiio  ililfereiitial  eipialion  for  7'.(?)  il  is  foiinii  that 
/.ero  heal  transfer  aloii<;  the  wall  can  he  ohtaineil  for  ii  ni>n-  | 
uniform  (as  wi'll  as  a  uniform)  wall  lem|)erature  ilistrihiilioii. 
'I’his  result,  ill  fact,  Generalizes  a  result  tif  refereiiee  47  tef. 
ref.  (>  for  iletails). 

P'rom  a  (iraelieal  (toini  of  view,  it  shoiilil  he  ki'jit  in  iniiiil 
(hat  the  solutions  ilevelo|)eil  here  are  haseil  on  the  viseosily- 
lemjieraliire  ri'lations  (ti)  anil  (7).  which  are  an  a|)|)roxima- 
tion  to  the  actual  relation  for  air.  Because  of  ri'laliims  tii)  : 
anil  (7).  the  results  ohtaineil  here,  namely.  ei|ualii>ns  (.’{o) 


anil  (4()a).  iiidieale  lhat.  for  a  fixed  wjill  I  einpeia  I  me  ih.’ 
skin-friel  ion  eoetiieieiil  and  (he  Niissell  numhei  will  he 
inile|)emlenl  of  Mach  nuiuher.  l'’tir  (he  Sullierland  \  i-- 
eosily-lem|)eral  me  relation,  however,  ihi'  will  nm  he  ipiiie 
valiil  (ef.  ref.  4S). 


KM)W  WITH  f'RKSSl'RK  (tRADlK.NT,  PRANDTI.  MMHKR  l‘r  I. 
AND  TNIKORM  WAI.I.  TKMI'KRATl  RK  >  li  t  ONSTAM 


P'rom  the  etpiations  derived  in  (he  section  "Bii'i'-  P.ipin- 
lions."  a  relatively  simjile  ami  siillieieiil l\  aeeinale  e  ihiid 
for  most  prael  ieal  |)ur| loses  of  ealeiilal  ine  1.1111111:1 1  hiimid.i  1  \  - 
layer  fharaeterislies  in  a  |)ressure  eradieiil  w  ith  heal  iraii'fei 
w  !1  he  tlevelo|)eil .  P'or  this  |)ur(>ose  il  will  he  a'Siinied  ihal 
the  Brandtl  niimher  of  ihe  lliiid  is  miily  and  llial  ihe  wall 
lem|)eral  lire  is  uniform.  'Phese  resi  riel  ions  eoii'idei  ;i hi \ 
sim|)lif\  the  nialhemalieal  :inal\sis  (ef..  e.  e..  e(|'.  i2I  '  ami 
(2211. 

.V  furl  her  advantiiGe  of  assiimiiiir I  hereisihal  iiilliis 
east'  il  ftillows  from  the  eiierey  |):irlial  dill'ereiil  ial  eipialion 
(.■>1  that  for  /.ero  heal  lr:insfer  ii.  e..  for  l07  di/e.  K  ami 
lienee,  (d/I  di/t.,  0)  II  (’olislanl  leGartlles'  ol  llie  prc'siire 
Gratlienl.  'Phis  \  ields  I  he  followiiiG  w  ell-know  n  value  ol 
the  et|iiilihriiim  adiahat ie  wall  tem|)eraliire  7',  lor  a  I’randll 
nuiuher  of  I : 

- '/’i  (  1 '  .1/,-')  /■  (  1  +  '  ,\/  )  '(li, 


(’onsetpielilly. 


7; 

7; 


i47 


'I'liiis.  the  parameler  A  is  in  this  ease  ihe  ph \  sie;ill,\  siGnitieaiii 
ratio  of  (he  aeliial  wall  leinperaliire  10  ihe  eipiilihriiim 
atliahalie  wall  lem|)eral  lire.  It  follows  from  this  iliai  for 
zero  heat  Iransfer  A  -l.  Il  will  he  seen  (ef.  ei|S,  I4.-1I. 
(.■).-)a).  ami  (.■).')hii  that  this  eomlition  is  exaelly  saiislieil  h\ 
the  a|)|>roximale  etpiations  and  solutions  used  here. 

.V  brief  discussion  of  methods  for  eases  of  /'/ 7^  1  and  or 
nonuniform  wall  tem|)eratiire  has  heen  ,Gi\en  in  the  intro- 
tliietitm.  'Phe  develo|)metit  ynveii  here  will  he  essentially 
the  same  as  lhat  in  references  ti  and  7. 


tiKNKHAl.  Al'I'RtlXIMATK  St>l.nTt)N 

With  /'/  I  ami  A  eonstalit,  while  i/,  n  Iji  is  :irhil  r:ir\ . 
the  etiellieieiits  11  „  and  A,,  in  etpiations  ( I.':.'  and  ( Itih  h\  \  irliie 
t)f  hoiimlaiw  eomlitions  (I!))  to  {2.tl.  can  all  he  ex()ressed  in 
terms  of  n-.  and  A,,  where  A,  remains  arhitr;ir\,  while  </.  is 
•riven  hy  : 


-Kt  ii' 


(4Si 


'Phe  (irtililes  in  terms  tif  tij  anti  A,  are  then: 


;(2r  or'  •  hr  '  2r")  ■  «».  ■■))|  ' 'Jr  It)-'  -  lOr  -•!/  • 

"1 

(A|t)A)(  -r-  lOr'  -2()r'e  lor’  '4t'')|  I41h 


II 

//| 


A  :  (I-  A)h>or'  S4r'  -  7(lr''-  2(lr'l  - 
A|(r  2()r'  ■  4.")r  '  -iftir’’  •  lOr') 


CiO) 


10 


12  1'. 


\.\  llitXAl,  Ain  ISdH^  CIIMMI 11  i.i.  liil;  Al  K(i\ At  I  l(> 


^Vitll  prolilcs  (4!l)  iiiid  tlic  folinu  iiiir  cxplicil  i-xprcssidiis 

fof  /'i,  and  /•,.  ail'  ()l)lain('d: 

F,  O.lOiia  •  O.DD'illa..  . 

(1.001)4  I2(  A, a.,  hi  ().00000i).‘)(/),a,  /()- 

o.oooir);{(/>|(/,-  hi 

/■’,  (). .40 ■)-  (), nOOl  1  /,)  •  0.1076,  • 

0.02r-*<(,  -  ().000(12(/..-  ■  0.002.st6,((.  Ii\- 

0.000l,')(6i(;,-’  hi-  0.000000.4(6, a,  6l-' 

F,  (I -/()|0.24)1  -  O.Ol.w,  O.OOl.Sl  (6|a,  6)]  - 

6|[0.0)).S4  0.00424a,  0.00041  (6|a,  6)| 

4Ailli  cxpii's.'^ions  (401.  (.40i,  and  (.41)  inscrii'd  in  lapialions 
(0)  and  (10).  two  oi'dinaiy  diH'cii'nl ial  ('(luations  I'oi  \(f)  and 
6,(f)  arc  olilaincd.  .\ltliou<rli  these  can  hesolveil  iinnierieall v 
for  a  ^iveii  dislrilmlion  of  a,  >i  Ai).  tlie  process  may  l>c 
tedious.  A  relativelx  sim|)le  jieneral  approximate  solution 
of  tliesi'  ('(pialions  will,  t lii'refori',  he  derived. 

K(|iiation  (0)  can  lie  solved  approximately  foi  X  In  assiim- 
inir  tlial  Fi  and  /•',  I'an  lie  I'eplai-ed  theix'  1)\  constant  “aver- 
aire”  values  and  /■',  over  the  distance  J.  This  is  Justified 
In  the  fact  that  the  variahle  terms  there,  nhieh  ari'  propor¬ 
tional  to  a,  and  6,.  are  relatively  small  (ef,  eqs.  (.‘ilii.  This 
is  equivalent  to  rephieinir  a,  and  6,  hy  eonstaiit  avenive 
values  (7,  and  6,  for  this  pitrpose.  With  (‘<|uation  (40)  for  tlii' 
\eloeity  |)rolile  and  equation  (4.S)  for  a,,  equation  (0)  then 
heeomes  the  follow  ilur  linejir  ordinaiw  dilferential  eipiation 
in  X: 

F;  2)X'  ■  X  /-'ll  Pi'  Pi  I  '  la,'  a,)j^v'i+  '  '  .1/,-  tv,  F,) 

27 '(p,  pi)  (  7',  7’, )  (a  .  a,)  (.">2) 

where  vi  is  a  eotistant  ydveii  hy: 

Vi  0.4/i  ■  t),t)00.46,  -0,tltl44.s  0. 0242a,  O.Ot)l24<7,- 

(0.t),s.4,s--0.(i(i4.'>,S(7,i(,7,6,  406)  (.->4) 

With  relations  (14)  and  (I4i.  the  solution  of  etpiation  (.421 
satisfyiiiyMhe  condition  X  t)  or  a  finite  value  (if  a,  0  at 
i  0)  al  the  leadini^  edire  ^  (I  is  found  to  he: 


\  ■reiieral  approxiimi I e  snlulioii  hu-  li,tA  can  he  uhlauud 
in  a  eompaiativel\  simple  manner  h\  l■linuualmy  X'  fiom 
dillenail  ial  e<|uatious  (Oi  ami  (Itli  in  eoujmwliou  with  the 
saiiK'  l,\pe  ol  simplilyiui;  approximation^  eiuii  ermii^  the  a 
terms  as  made  in  ileriviu”;  equation  1.44'.  i.'ms'  ii-i'.  q  hu 
details.)  A  (piadi'alie  e(|ualion  in  6  is  therein  Khlaiiud. 
with  the  siilulion 


'l, 

—  ./ 

'(./  ) 

/ 

'T 

'J/// 

'  J//I  / 

7/6 

w  here 


III  ~A  10  '  )i4.t)()  l.  P.ttia,  4  42.70.4  0 


2llt)0a  1 


/■  0.24002  ■  a. It). 070!  7  O.OO.I.so.v,; .  i  y  |o 
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^  2.S4..' 


1  ■  (42.20  2.0(i7a  la,  i7.000  (I.Od.sO. 


/  -2(1  6|((l. 24002  0. 01400a 


0.46  :  a, (0.021  10  0,0000: 


la,)  I  !+';;[ 

I02l0a,ijy 


0.  IOfO.4 


The  phvsieally  a ppi’opritl I e  root  in  equiilion  (.4,4!ii  will,  m 
ireneral.  la'  lluil  which  is  closer  to  the  value  2(1  In. 

After  Xf^)  has  heeti  ohtaiiied  h\  means  uf  e(|uation  .41  . 
the  eoedieieui  a,(i)  follows  fnm)  efjualion  (4Si.  and  64-.  e.-oi 
tin'll  he  direell\-  ealeilhileil  hy  means  of  equations  (.4.4i, 
For  ohjeels  with  sharp  leadiiev  edjres,  for  w  hich  X  0  al  i  0. 
it  will  orilinai'ily  he  found  that  an  appi’oximale  value  of  6 
aeeordinjr  to  equations  (4.41  is  lhal  eiveii  h\  equation  !4.4., 
which  is  valid  exactly  for  the  ease  a,:=().  This  is  ilhisl  raleil 
in  detail  in  reference  7  hy  numerical  exanqile  I'm-  the  supi'r- 
sonie  How  over  a  thin  hieonvex  airfoil. 

The  aeiieral  ap]>ro.\imait'  sohitions  oiveii  h\-  e(|uaiions  i44- 
and  (.4.4)  are  (piile  eonveiiieiit  for  aeliial  ealeiilalions  and 
invohe.  al  most,  numerical  int(';j:iation.  These  solutions 
will  he  approximately  valid  as  lone  as  the  11.  terms  in  expiis- 
sions  (41)  ioe  indeed  relatively  small  either  indi\iiluall\  or 
collect ively.  Such  is  expected  to  he  ordinariK  the  ease  in 
practice.  In  eases  for  which  the  a,  terms  heeonie  rehilixi'h 


I  ('/i  II , )  '1  (7’,  7' A  '  'o/s 

X  i4  6,)(  ‘  "  ,  -,.,1  (.44) 

'"1  11 J  ''  t'/Y  Tj'  ‘  'f' 

Kipiation  (.44)  is  simiho'  in  form  to  e(piatioiis  oliiaineil  for 
zero  heal  transfer  in  references  2  and  4!)  and  to  lliose  oh- 
tailied  for  heat  transfer,  hill  with  fourth-decree  profiles  and 
twd  homidary-laver  thieknesses.  in  reh  reiiees  I!)  to  21.  It  is 
intereslinc  to  ohserve  that  for  zero  heal  transfer  it  is  possible 
to  derive  forms  like  eipiation  (44)  (ef.  refs.  2fl  and  451)  hv 
appliiti};  the  .Siewarf.son-Illincworth  transformation  (refs. 
.40  and  41 ).  Ilowi-ver.  h\-  tin-  [ire-senl  (a|)pro,\imate)  mefliod 
of  analysis,  it  is  seen  that  w  ith  the  nseof  only  the  Dorodiiitsyn 
transformation  (.S)  such  a  form  can  he  straightforwardly 
derivi'd  even  for  the  ease  of  heat  transfer,  lint  nniforin 
temperature,  alone  the  wall  (ef.  also  ref.  20). 


larce,  however,  the  ordinary  dilferential  equations  (lii  and 
III))  ma\  have  to  he  solved  numeiiealK  . 

In  exalnatini;  /',  and  vi  a  reasoiiahle  avera<;e  value  7  for 
;  II:  for  an \-  civen  n,  ii^(^i.  6,  and  .1/.,  can  usually  heohiained 
,  hy  eonsideriiic  eipiation  (4.S)  for  (ii,f6X)  and  equation  (.44! 

I  for  \  ('.  A  satisfactory  averayu'  value  6,  for  6,  in  evahialiuc 
I  F,  and  vi  will  ordinarily  he  that  civen  hy  equation  (44). 
j  In  reference  2.  mimerieal  examples  hased  on  the  ease 
i  II, 'll:,:  -  I  (where  6  is  a  positive  eonslant)  for  .16,  0.  I. 

j  and  4  and  zero  heat  transfer  at  the  wall  (6  1 )  were  carried 

!  out  to  determine  the  aeeuraei'  of  ap|)roximate  solution  (.44) 

I  of  ordinaiw  dilferential  equation  (42).  (’oinjiarison  of 
I  the  .solutions  ohiained  hy  means  of  equation  (44)  was  made 
j  with  nnmerieal  solutions  of  dilferential  eipiation  (42)  with¬ 
out  the  use  of  any  of  the  appro.xiinatinc  assumptions  made 
in  derivinc  equation  (44).  'I'lie  eomparison  indicated,  on 
the  whole,  .satisfactory  acfeement  for  practical  purposes 
(ineluilinc  slahilit\-  ealeiilatious)  hetween  the  results  of 


ANALYSIS  ol  LAMlNAli  fOMl'KKSSIBLK  Ht)l  \l)AHV  I.AVKK  WITH  HKAl  IKANSl  l  l!  AM)  Fl!l»l  liK  (.liADlKV  I 


(<(|iiiili(>n  (54)  1111(1  (lie  niiiiK'l'ioil  solution  of  (‘(jiiiil ion  (.yji.  ; 
Dctiiils  )tn'  <;iv(‘n  in  reference  2.  Similar  eomparisons  : 
liiive  also  l)e('n  earrieil  out  in  refeiA'iiec'  21  for  tlie  eases  i/,  a  ,  ] 

1  :L:/)t\vith  heat  tratisferat  the  wall.  The  ajrreenieni  hetween  ! 
the  t  \  pe  of  approxitnate  soltttion  fjivt'ti  h\  e(|tiatioti  (■>4)  and 
the  tnitnerieal  .sohitioti  of  the  origitial  orditiarv  difIVrenlial 
('(ptation  was,  apiin,  foiitid  to  he  oti  the  whoU'  satisfaetorv.' 

SKIN  PHICTION.  HE.tT  TR.4NSKEK.  VELOCITY,  .YNI)  TEMPEHATI'HE 
PROFILES 

With  X({)  and  A|(^)  determined,  the  hoimdarv-la\ ei- 
eharaeti'i'isties  can  all  he  straijilitforwardly  ealeiilated.  The 
local  skin-friction  eoi'llieient  will  he 

_(m  dll  di/l„ 

"""  1 
2  ^  "  ' 

4|i  (It:  .■)!-  (/),(/j  (r  7', )()),  1/ , (.Ati) 

The  .\nsselt  mimhi'f  ^ivinjj:  local  heat-trjinsfer  properties 
iit  t he  wall  will  lx- 

The  \  ('loeit  v  and  temperature  profiles  follow  from  e(|iiations 
il'.il,  and  i24)  in  eonjnnetion  witli  e((tiation  IS)  for 

transforming  to  tlu'  physiciil  plane.  Kor  /.ero  heat  liitnsfer 
at  till'  Willi,  an  ex|)lieit  exprc'ssion  for  i/  as  it  fimetion  of 
r  in  terms  of  n,  is  eiven  in  iippeiulix  .V  of  reference  2.  'Phis 
expression  can  lie  eonveiiientlv  written  in  the  form 

'/It.  r)  ,  /  7--  1  \  .  t 

-  .)  ) -'A'  U'l  (  r)  •  r)  •  T)| 

where  '.'.i,  itnd  gj  iire  definiti'  fimetions  ( |)olynomiitls)  of 
r  oidy  which  remain  the  sitme  for  itll  eases.  These  fimetions 
can.  if  desired,  he  evaluated  itnd  plotted  once  for  till.  A 
similiir  expression  can  lie  obtained  for  the  ciise  of  hciit 
transfer  itt  the  witll,  except  thiil  iidditioiiiil  terms,  such 
iis  tliose  proportional  to  (/j6|  h.  will  he  included.  'Pile  imi- 
versid  fimetions  of  r  thus  obtained  eiin,  if  (h'sired.  also  he 
eviiluiited  once  for  all.  For  it  eiven  vahi(>  of  y  or  y  5 
tef.  e(|s,  (21)11)  and  (2t)h)  )  can  then,  in  any  ^ivett  ease,  he  i 
found  (piite  straiglitforwiirdly  for  values  of  r  from  r  -0 
tori.  i 

,V  numerical  example  to  cheek  the  aeettrae,\  of  the  results  ^ 
obtained  by  tlu'  ecpiations  di'veloped  in  this  section  was  | 
carried  out  in  reference  2,  'Phis  example,  as  previously  | 
indicated,  was  the  ease  of  (low  with  a  linearly  deereasiiifr  : 
velocity  outside  of  the  homidarv  layer.  X'elocify  profiles,  1 
local  skin-friction  coellicient,  and  minimum  critical  Reynolds  : 
numbers  for  laminar  instability  were  ealeiilated  by  this  ' 
means,  and  the  results  for  incompressible  flow  (/)  “!  and  j 
0)  thus  obtained  were  compared  with  those  based  on  j 
the  series  solution  in  refereiiee  r)2  of  the  orifrinal  partial 
difl'ereiitial  equation  (I).  'Phe  agreement  was  in  all  cases 
found  to  be  satisfaetorv  for  practical  purposes.  (Details 
are  given  in  ref.  2.) 

'Phe  solutions  presented  here  require  some  modificatiou 

»  In  Tfl'Tt'nrp  '^l,  thi-  unrtlysis  wiis  bnvi’il  on  profih-s  in  oonjtiftrlion  with  !M»th  a 

ffjiTtmil  f»n*l  a  lynainicHl  hoiimlary-laycr  thUiknt'.ss.  An  ctiijation  tinite  analoeons  (o  rtiua- 
lir)n  f.'it'  for  ^f^),  howt'ViT,  whs  obtiilnptl  by  finalok't)iis  approxlrnatin;:  tiss'iniptions. 


in  two  important  special  cases;  lai  Flow  near  a  forward 
stagnation  point  and  (hi  calcnhilion  of  the  separation  point 
in  an  adverse  pressure  gradient 

STAONATlON  EI.OWs 

Phe  case 

a,  II  /)s  .‘is 

where  li  is  a  positixe  eonslani  represents  phx'icalK  the  How 
in  the  vieinitx  of  a  forward  stagnation  point,  siah  a~  the 
siibsonie  (low  over  the  leading  edge  of  a  blunt  object  For 
zero  .\laeh  number,  an  exact  solution  of  the  ordinarx  dilfcr- 
eiitial  eipnitions  (tl)  and  lit))  ixvith  uniform  xx:dl  lempcra- 
tiire)  can  be  obtaiiu'd  in  tlie  form  X  ('onsiani  and  /< 
Coiislant.  Kcpiations  (!))  and  (lll>  then  become  algebraic 
eipiations  for  X  and  li,.  For  the  special  case  of  zero  hc;ii 
transfer  t//  1  and  li,  11).  howcxi'r,  it  has  alreadx  been 
found  (refs.  4d  and  2)  that  these  eipiations  w  ill  not  x  leld 
any  physically  signilicant  real  roots.  In  reference  2  it  wa' 
slioxvn  that  an  approximate  solution  can  still  he  obtained 
ill  this  case  by  xvriting  the  algebraic  eipiation  as  /iX'  h 
and  taking  the  value  of  X  for  xvhich  /( X 1  has  a  local  'naxiniimi 
value  relatively  close  to  the  X-axis,  'Phe  root  X  h.lsl  wa- 
thus  obtained.  'Phis  solution,  hoxvex’er,  is  iinsal  i'faclorx 
in  principle.  ( 'onse(|uenlly ,  the  use  of  fourth-degree,  instead 
of  sixth-degree,  xelocitx  profiles  for  this  ease  was  inxesli- 
gated  in  refereiiee  2.  'Phe  proliles  xvere  chosen  to  satisfx 
the  usual  Karnian-l’ohihausen  conditions,  \  phxsicallx 
signilicant  real  root,  namely.  X  7.il.'>2.  was  now  obtaineil, 
and  the  aeeiiraex’  of  the  resulting  solution  xvas  compared 
xvitli  till'  results  of  an  exael  solution  iref.  .all.  In  particular, 
skill  friction,  velocity  profiles,  and  minimum  eritieal  Keynolds 
nnmber  xvere  compared.  'Phe  eomparison  indicated  that  the 
results  obtained  by  the  use  of  the  fourth-degree  proliles  led 
on  the  xvhole  to  results  of  satisfactory  accuracy.  It  xvas 
therefore  eoneluded  that  the  bomiihirx-layer  charaeteristics 
in  lloxv  near  a  forxvard  stagnation  point  can  be  ilelermined 
xvith  satisfactory  aeciiraey  liy  the  Karman-l’ohihaiisen 
iiiethod  xvith  fourth-degree  proliles. 

'Po  ealeiilale  the  boundary  layer  near  a  forxxard  siagnaiion 
point  for  the  more  general  ease  of  heal  transfer  at  the  wall, 
in  particular  for  a  preseribed  uniform  xx all-teiiiperat me  ratio 
/)  or  1\^,  the  method  ol  reference  2  can  be  generalized  bx 
inlrodiieing  fourth-degree  stagnal ion-eiit halpy .  as  ell  as 
velocity,  profiles.  'Phis  has  been  carried  out  in  reference  20 
wit  h  the  ini  rod  net  ion  of  a  t  hernial,  in  add  it  ion  to  a  dynamical, 
boundary-layer  lliiekiiess.  'Pxvo  algebraic  equations  in  essen¬ 
tially  the  Ixvo  (coiistiint )  boundary-layer  thicknesses  are  ob¬ 
tained.  'Phese  e(|uati()ns  can.  in  general,  be  solved  either 
iiuiuerieally  for  a  given  h  or  by  using  the  values  in  figures  1 
and  2  of  reference  20.  Although  reference  20  is  based  on 
lloxv  over  a  sxveat-eooled  surface  and.  lienee,  ineliides  a  nor¬ 
mal  mass  lloxv  at  the  wall  (r^  r„  at  r  0),  the  results  there 
can  also  be  used  for  an  imiiermeable  xvall  bx  simply  imlliiig 
f'=0  and  letting  h  be  arbitrary.  ('Phis  ('  is  not  to  be  con¬ 
fused  with  the  temiieratiire-viseositx  factor  used  in  the  ])res- 
eiit  paper.)  An  e.xam|)le  of  loxv-siieed  (.1/^=0)  lloxv  in  a 
fax’orable  pressure  gradient  xvith  a  stagnation  point  at 
J  -0.  representing  siibsonie  lloxv  over  a  turbine  blade,  xvas 
carried  out  in  reference  20  on  the  basis  of  the  method  pre¬ 
sented  there. 
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UKHOKT  r24:>  —  NAl'IONAL  AUVISDHV  COMMl  11  Kl-,  lOU  Aim  IN  Al  IK'S 


All  alti'i'iiiilivc  iiu'lliod  of  calciiliU iiijr  Hows  iioar  a  forward 
slaoiiatioii  point,  liascal  on  llif  use  of  a  sinolo  Ixnindarv-laycr 
lliickncss,  is  oivcii  in  roffrcnco  22. 

<  AlA  in.ATION  or  -SKPABATION  POINT 


'rill'  I'linations  tlins  far  developed  in  tins  section  can  he  , 
used  to  ealeillate  the  laminar  sejiaratioti  point  in  an  adverse 
pressure  >rradient.  The  results  thus  obtained  will  frenei-allv 
he  tnore  aeenrate  than  those  obtained  hv  the  use  of  fourth- 
decree  proliles.  By  an  analysis  for  ineoinpri'ssihle  How  foi' 
till'  ease  of  a  lini'arly  diininishing  velocity  outside  the  hound-  j 
ary  layer,  it  was  fonnd  (ref.  4.")).  however,  that  still  <ri<'ater  | 
aeeitraey  for  the  location  of  the  separation  point  is  obtainable  ■ 
h\'  the  special  use,  for  this  [iiirpose,  of  .seventh-decree  velocity 
profiles  satisfvinc  an  additional  condition  involvinc  the 
fourth  derivative  of  the  velocity  at  the  wall  at  the  separation 
point.  This  condition  wonld  necessarily  hi'  satisfied  by  an  i 
e.xaet  solntion  of  the  oricinal  partial  differential  ec|nations. 
'Phis  method  of  ealeulatinc  tlii'  separation  point  was  nuIi-  ' 
seipienth  I'xtended  to  eoiupressible  flow  with  zero  heat 
transfer  in  reference  2  and  to  compressible  flow  with  heat  , 
transfer  in  reference  7.  The  method  of  analysis  to  be  pri'-  i 
si'iited  here  is  essetitially  that  of  referenei'  7. 

It  tmiy  be  metitioneil  that  a  considerable  mimbi'r  of  ! 
methods  of  ealeulatinc  the  laminar  separation  point  have  : 
been  developed.  Xo  attetnpt  will  be  made  here  to  summarize 
or  evalnati'  all  of  tiu'si'  methods.  For  incotnpri'ssible  flow, 
a  mi'thod  which  has  heen  foitnd  to  yield  results  of  satisfae-  ■ 
tory  aceitracv  in  addition  to  that  of  reference  4")  is  that  of 
reference  .■)4.  For  cotnpressible  flow  with  zero  heat  transfer  ■ 
(which,  of  course,  includes  incotiipressible  flow)  recent  nu'th-  : 
oils,  in  addition  to  that  of  reference  2.  are  those  of  referi'tices 
20.  ."lO.  ■).").  and  .7(),  For  cotnpressible  flow  with  heat  transfer, 
ihe  only  studies  of  lamimir  separation  which  apjiear  to  have  ; 
been  tiiade.  in  addition  to  that  of  referenei'  7,  are  those  of  i 
references  b!  to  It'i.  'Phe  advantace  of  the  method  to  be  I 
presetiti'd  here  is.  once  acain.  not  onl\  that  it  appears  to  . 
yield  results  of  adi'nuate  accuracy  but  that  the  analysis  is 
kept  relatively  sim|)le,  althouch  it  is  based  on  a  minimum  of  ; 
what  micht  be  termed  mathematically  “arbitrary”  assump-  j 
tions.  The  method  of  analysis  developed  heri'  is  indeed  | 
suffieieiitly  simjHe  and  fle.xible  to  be  applicable  to  a  wide  i 
variety  of  eonditions.  ('Phe  method  has,  in  fact.  Iieen  (piite  i 
recently  extendi'd  to  thi'  case  of  eom|)ressible  flow  over  a 
transpiration-cooled  surface  (ref.  .'iP).)  'Phe  calculations  to  i 
be  performed  aeeordinc  to  the  method  presi'iiti'd  here  will 
be  relatively  simph'  and  will  involve,  at  tnost.  numerical  j 
inti'cration.  ! 

By  difrerentiatinc  the  mometitum  partial  differi'litial  ! 
eijuation  (1)  it  eati  be  shown  (ref.  7).  undi'f  the  pre.sent  j 
assutn|)tion  of  a  Brandi  number  of  1  attd  a  litiear  viscosity- 
tem|)erature  I'l'latioti.  that,  at  the  separation  point,  with  or 
without  lii'itt  tratisfi'i- at  the  wall. 

(db//dB)„=0  (.Tft) 

'Phe  sevi'tith-decree  veloi'ity  jirofile  satisfvinc  condition 
(5!))  in  addition  to  cotiditions  (17)  to  (23)  is  “ 

*  'rhb  st‘VBnfh-«li*Rr(‘o  iiritfllB  Is  to  ii.st'tl  only  for  polculntton  of  the  sf’jmration  For 

other  (nirpose.s,  the  slxth-'lei»n'e  profile  fo<j.  (IIP)  shotil‘1  be  tirotl,  even  In  an  ivlverse  itrepsure 
'.'ra»IUTt  <cf.  ref.  2i. 


where  <ij  is  ci'eii  by  equation  i4Si.  Separation  ocem- 
where  {dn  d;/),,  I)  and.  hence,  where  {  dn  Or  i,,  1).  'Phcrc- 

fore.  aeeordinc  to  equation  IliDl,  separation  will  oenn  ulicii 
has  the  value  (denoted  b\  ij,.; 


/|  +  .  - 
1 .) 


nl 


From  equations  (tsianil  (lill  it  follows  that  the  value  A 
id'Xat  the  separation  point  will,  in  cencral,  be 


A  satisfactory  approxitnation  for  /i,  in  e((ualion  lii2)  will, 
in  ci'iieral.  be  that  ci'en  by  equation  (4.'))  (cf.  also  ref  7' 
With  this  expression  for  /q.  eiiuation  (li2)  becomes 


X  10-r  1 

,71  ,,mll/i  - 4 


("i'".)(l-f  .,  ■l/i-) 


By  insertinc  prolile  KiO)  into  ililfereiitial  equation  lUi  and 
assuminc.  in  llu'  fori'coinc  analvsis.  that  the  and  h 
terms  in  /■',  and  F>  may  be  replaced  bv  constant  values,  an 
ordinary  ililferential  equation  of  the  same  form  as  eipiation 
(.j2)  is  obtained,  exi'ept  that  the  explicit  expressions  for  J'\ 
and  v?i  (to  he  denoted  now  as  .  and  ^i,.)  are  modiiied.  w  hile 
the  factor  2  on  the  richt  side  of  equation  (.■|2)  is  replaced  hv 
7  4.  ('omparison,  aeeordincl.v  .  with  the  solntion  (eq.  (.'ill' 
of  equation  (.■)21  yields  the  followinc  solution  for  Aiji  (de¬ 
noted  now  as  X, ) : 


-  r*' 

(fi  -- 


'Pakinc  the  constant  value  of  n,  as  that  at  the  se|)aration 
point  (as  in  refs.  2  and  7).  the  espressions  for  A',,  and  v”, ,  are 
found  to  he 


A|,=^().l  I  ().0()2o2.')a,,  I). OtI  1454(1... - 

(l.()()()0572(6,(/o,  /i)- -0.000574 (6|(r,,-  h  )  ■ 

O.OOOS.S7(6,(t..,7i) 

f  0)5) 

sPi..s-().2.5/H  0.0437-|-0.073.S/>,K).t):i4.S(t,,-- 

0.002t)l(fj,M-0.0077;i(6i"-'<  /')-0.001 147(/),aj;-'  A) 

0.0001 145fA,(/,,'/i)- 


where  Hj,  and  A,  are  ci'''’ti  by  eipiations  (01)  attd  (45).  re¬ 
spectively.  'Phe  qttantities  A',  <  and  .^i.  are  functions  of  A  oidy 
and  are  shown  in  ficure  1. 


AWI.ISIS  OK  LAMlNAIi  COMPHKSslHKK  HOINDAHV  LAMJli  Wllll  H  K  A  I  riiAN>lLlt  AM)  lMil»ri;l,  i.liADIIVI 


I  A 


For  any  jiivcn  rcfcronco  Macli  nuinlx-r  and  iiniforin 
wall  ti-mporaturc  ratio  fi.  tlio  separation  point  in  a  re<:ion  of 
iriven  adverse  pressure'  <;radient,  as  sjeecified  ley  u, 
will  l)e  tlie  station  ?  at  evliicli  tlie  rijrlit  sides  of  eqiiafiotis 
{^y^)  and  ((>4)  are'  eepial.  'I'lnis.  it  is  nee'e'ssary .  in  fte'iie'i'al. 
otdy  tee  jeleit  X  ve'rsus  in  the'  antie'ipate'el  vie'iiiity  e>f  se'para- 
tion,  in  ae'e'eerelane'e'  with  both  eepiatieuis  Hi'A)  anel  (fU)  anel  te> 
eh'te'riiiine'  the*  pe>int  of  iiiterse'ction  of  tlie'se'  two  curve's. 
'I'he  se'paration  [)oint  Avill  evielently  he'  iuele'|>e‘ueie'nt  eef  se> 
that  for  the'  |)urpose  eef  eleteriuiniusr  the'  se'paratieui  point 
one'  may  se-t  t'=il. 

In  case  the'  re'frieui  e)f  aelve'ise'  ()re'ssure'  jrraelie'iit  starts  at 
senile  peiint  elownstreaiii  of  the'  leaelinj:  e'il<re,  e'<(uation 

(t)4)  e-an  still  he'  applienl  elire'e'tly  in  I'ah'ulatinjr  the'  separation 
point.  (Ireate'i-  ae'e'urae'y,  heiwevi'r.  niifiht  he'  olitaiiu'el  in 
sue'h  a  e'ase'  hv  applying  i'(|uation  ((>4)  einly  feir  the'  re'gient  eif 
aelve'ise  pre'ssure'  graelienl.  Feir  this  purpeise',  e'i|uation  (t)4) 
must  he  moelifie'il  to  satisfy  the-  hounihiry  e'onelition  X  =  X,  at 
Thus, 


2->  - 1  •r,. 


wlie're 


^  >+i  '^1, 


(titi) 


and  whiTi'  X„  e-an  he  ohtaine'il  as  the  value  eif  X  at  ^  =  based 
on  <‘<(uation  (o4)  for  the  ri'gion  of  favorable  pre'ssiire 

gradient. 

For  purpose's  of  I'uleulating  the  separation  point  for 
various  value's  eif  the  temperature  ratio  7'„/7’«  and  of  the 
reference'  Mach  nunihe'r  er|uation  (dH)  may  lie  replaceil 


hy  'he  f))!l,)wiiig  ee|uivali-nl  e(|uali(in  for  llw  ealiic  of  X  al 
the'  separation  point  : 


lOof 


|4  - 11(7:, 


I  "  ■ 

'I' 


Fe|uatie>n  (1)7)  feilleiws  from  e)|Ualion  (eiM)  hy  in^erlio”  i.  la- 
tions  (14),  (4l)).  anel  (47i  llcre'.  with  1.4. 

.Nunie'iie  .el  e'Xiim|)le's  for  lle>\\  with  a  linearle  elci  i  c-.i'iii); 
ve'loe'ilv  at  the'  emte'r  e'elge'  of  t  he'  henindaiA  layi'f  are  ill  list  reil  eil 
in  eletail  in  re'fere'ne-e'  7.  anel  ihe'se  will  he  discussed  hrie  lh  iii 
the'  folleiwing  si'ction.  .\n  e'Xiiinple  hjisid  on  a  siagnalioii 
flow  followeel  hy  an  aelve'i'se  pre'Ssure'  giiidie'iil  is  also  eliscussed 
in  eh'tail  in  re'fere'iice'  7.  Foi-  ihe'  case'  eif  a  hnearl\  dccrciisin^ 
ve'leie  ity  ontsieh'  of  the'  heniiulary  laye'r  with  zi'iii  liciit  tiaiis- 
fe'r  at  the'  wiill  the'  se'|)aratioti  peiiiil  wees  i-alcula ' '  ■ !  lie  tic 
nietheiel  pre'se'tite'el  he'te'  for  .Mach  nuinhers  .1/  from 

II  le>  10.  The'se'  r'e'sults  are'  conipari'd  in  table  I  -c  of 

refe're'iie'e'  .‘lO,  atiel  the'  agreenii'tit  is  si'i'ii  to  inch 

I'h'se'. 

It  may  he'  ri'cidh'd  thill  the’  ini'lhoel  of  I'ahuhit iiiy  llu' 
se'paratietn  point  pre'si'iiti'el  lu're'  is  hiiseel  on  llii'  iissiiinpt  ion 
e)f  a  lini'iir  viscosily-te'inpe'iiiluri'  relation  moc /'  ainl  of  a 
Pranelll  nnmhe'r  /V  of  nnily.  It  is  nole'w  enl  he .  in  this 
e'onne'e't iem,  that  it  has  he'i'u  I'oni'lneh'el  in  a  rl'l■l‘nl  iinal 
(ref.  '>{))  that  for^ix7'*'  iiinl  ws,  1.  i he  si'pariit ion  point  <1 
/V>0.7  ex'i'iirs  a(  roughly  the'  same  position  iis  for  /'/  «,■  I. 


DISCI  S.SION  OK  SKIN-KKK  TION.  HEAT-TR.ANSKKK.  .SKP.AK.A- 
TION.  .AND  .STABILITY  (  H ARA(  TEKISTK  S| 

To  conclmJi'  lljj.s  repoj'l.  ii  suimmirA  will  he  gi\en  in  this 
se'e'tion  e)f  the  iinplieiit ions  of  the  I'ipnitions  ile'Se'leepeel  here 
'  re'gareling  the  elfeel  e)f  wall  leniperiilure,  Mach  nninher.  iind 
pre'ssure'  graelie'iit  eui  the'  laniinar-heninehiry-laye'r  l•haracle■r- 
istics.  The'se  I'oni'lusiems  have'  he'e'li  eh'ri\'e'il  anel  illustrated 
in  eh'lail  es|)eeially  in  re'IVrene'e'  7. 

SKIN-KKICTION  AND  H EAT-TRANSKKR  CelErKII'IEN'I'S 

The'  i'lfe'ct  of  wall  te'inpe'ral lire'  on  the  skin-friction  ami 
hent-tninsfe'f  i-oe'llicie'lits  will  eli'pi'liel  eui  the'  nature  (fiivorahlc 
e)r  aelve'i'se')  of  the'  pre'ssure' graelie'iit .  This  feilleiws  freun  the 
;  fae't  that  in  orelinary  elilfe'i'ent ial  e'epialion  (.■)2)  ami  in  expre's- 
sioii  (4S)  fell'  a,  the  te'mpe'ralure'  paranie'te'r  A  appe'ars  pri¬ 
marily  in  a  form  inulliiilie'el  hy  the'  ve'ieie'ity  graelie'iit  a,'. 
;  The'  effe'e't  elf  the'  wall  ti'iniii'ratiire'  em  the'  skin-fiict ion 
!  I'eie'flie'ie'nt  arising  freun  the'  (Wi'/i)  te'rm  in  e/j  is  part ie  ularly 
important .  Thus,  eepiations  (.■)())  ami  (4S)  sheiw  that, 
I  witiiout  the'  I'ffee't  eif  the  te'inpe'rat uri'-viscosit \  fiii'tor  < 
i  lowering  the'  wall  tempi'rature'  ti'iiels  to  eliniinish  the'  leii-al 
I  skin  friction  in  a  favoralile  |iri'ssuri'  graelieni  (negative'  a,') 
i  anel  to  increase  it  in  an  aelverse'  pre'ssure'  graelie'iit.  It  I'an 
he'  shown  (re'f.  7)  that  a  similar,  hut  inncli  smalh'r.  I'ffecl 
on  the  Nusselt  nuinlier  will  a  ii'iiel  to  oi'i'Ur. 

Since  the  velocity  graelient  a,'  in  the'  i'(|uatieins  eli've'lei|ii'il 
here  (cf.,  e'spe'e'ially,  I'ejS.  (4.S),  (.72),  anel  (.72)1  appe'ars  in  a 
form  multiplied  In'  the'  wall-le'inpi'raturi'  ratiei  h.  it  I'liii  he 
inferreel  that  a  lowe'l'ing  eif  the'  wall  ti'inpe'raturi'  has  a  te'iiel- 
ency  to  diminish  the'  elire'e't  elfe'e-t  eif  a  givi'ii  pre'ssure'  graelie'iit, 
that  is.  the  effi'et  of  a,'  as  such,  011  the  houiulary-layi'r  iirei|i- 
erties.  This  is  e'.xphiinahle  physie'iilly  hy  the'  ini're'ase'el 
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iiiipoitiincc  Ilf  llu‘  iiirrlia  foiri's  ii‘lati\i'  lo  llic  pii-ssuie 
irinilii'nl  Ix'CiUlsi'  of  llu'  iiicri'iisc  iif  llic  lluiil  ili'iisily  whan 
till'  Willi  ti'iupciiitilia  is  ilcciciisi'il.'"  A  i‘lciir  illiisi  rat  ion 
of  this  will  ha  saaii  siihsa(|Uaiil  l\  in  lha  iiiial  \  sis  of  hmiiiiar 
saparalioii  It  must  ha  ohsarvad.  howcvar.  that  lha  ali'aat 
of  ii  prassiira  iiiiiiliaiit  iilso  iippaiirs  iinlirai  tl\ ,  iiamal\ ,  in  the 
variation  of  a,  a,  iind  l\  T.  with  For  Maah  niimhars 
iihova  I.  in  fact,  tha  7’,  /',  tarins  in  \  {aip  f.")4)i  may  haaoma 
part ianlarly  impoitant.  so  tinit  in  snah  a  aasa  tha  nal  alli-al 
of  tha  prassiira  aradiant  may  aatnall\  ha  inaraasad  hy  a 
lowarina  of  tha  wiill  tamparattira.  'I'his  is  illiistriilad  in 
datiiil  in  rafaranaa  7  hy  ii  niiniaiiaal  axampla  for  lha  siipar- 
sonia  How  over  a  thin  tiirfoil  (aspaaially  at  .1/,  d). 

From  aipnUions  t.'ilii  and  toTl  it  follows  that  tha  ratio  of 
local  skin  friction  to  Xiissalt  nnmhai'  ciin  ha  axprassad  in 
tha  form 


A  //  V  ' 


F'or  How  iilonir  a  Ihit  plate  ( a,  a  ^  1 ,  a^  tl.  and /),  2(1  -//)). 

aipiation  (liSi  implies  that  /  -‘J.  Foi  How  in  a  prassiira  | 
aradiaiit.  howavar,  since  ordimirily  /i,  ~ 'if  1  — /< ),  it  follows 
from  aipnilions  |li,S)  and  (4S)  that  idoni;  lha  How  in  a 

favorahlc  prassiira  vradiani  (ai'>  l)).  and  /  <  2  in  an  adverse  . 
prassnra  ^radiant  (a,'<lli.  Moraovaf.  it  also  follows  from 
lhasa  aipiations  that  lowarinj;'  tha  w all-talitpaialufa  paramatar 
/(  will  tend  to  hrinv  a  alosar  to  its  value  for  How  without  a 
prassiira  irradiaiit.  This  illustratas  tha  diminution  of  tha 
diraat  ali'aat  of  a  prassiira  >riadi<‘nl  hy  aoolinir  of  the  wall. 

From  aiiiiations  (541,  (5lii.  and  (57)  it  follows  that  liot/i 
tha  skiti  friction  and  .Nitssalt  niimhar  will  ha  proportional  to 
\  t".  Thus,  an  aH'aat  of  wall  tamparattira  on  tha  skin-  . 
friction  ami  haat-traiisfar  coaHiaiaiits  follows  from  lha  vis- 
aosity-lam|)arat lira  aoaHiaiaiit  f  arising  from  tha  partiaiilar 
visaosilv-taniparatura  rahilion  (a<|s,  (til  and  (7))  assuiiiad 
hare.  This  ali'aat  is  indapandaiil  of  lha  prassiira  >rradiaiit. 
From  acpiation  (7)  it  follows  that  if,  as  will  ordimirily  ha  ; 
tha  aasa.  T„'';  S.  that  is,  7',, '7  211')°  K,  than  a  low ('riiiir  of  tha 
ratio  T„  l\  will  inaraasa  and  .Va.  For  a  fixed  ratio  ; 

Ii  of  /'„  /',.  it  follows  from  aipiation  (25)  that  a  Maidi  miinhar 
aH'aat  will  also  appear  in  Thus,  if  7'„(>.S’,  than  for  a 
fixed  value  of  h  an  inaraasa  of  Maah  nuitihar  M ^  w  ill  diminish 
('  and  liaiiaa  will  land,  as  far  as  ('  is  aoiiaariiad,  to  diminish 
hoth  tha  skin-frialion  and  haat-t ransfar  aoaHiaianIs  in  pro¬ 
port  ion  lo  ^  4 '.  ; 

From  acpiation  f54l.  as  has  alraad,\’  haali  liotad.  it  will  ha  ! 
foiinil  that  in  the  |)rasanaa  of  a  |)rassuii'  ^radiant  \  ('  may  ha 
appraaiahly  affaalad  hy  tha  Maah  numhar  haaaiisa  of  the-  ; 
values  of  T,  7\(^).  Consacpiaiitly,  it  aan  ha  infarrad  that  a  i 
prassiira  vradiani  will,  in  >r('naral.  land  to  anhanaa  tha  I'fl'i'al 
of  .Maah  numhar  on  hofh  lha  skin-friction  and  haal-lransfi‘r 
aoaffiaialtls.  This  affaat  will  da|)and  on  tha  natiira  of  lha 
prassiira  f;radianl.  For  a  favoniltia  prassiira  fjnidianl.  for 
axampla.  for  whiah  )/,'ii^>l  and  haiiac'  7’,,'7'„<(1,  an 
inaraa.sa  of  .Maah  niimhi'r  will  tend  lo  inaraasa  X/f  and, 

'I  hN  rttm  liisutii  fintl  thi-  ffiri  '.'oinu  cuiKliistdris  on  fht-  flTfcT  of  wnll  toniivrsitnn-  on  skin 
friction  have  liot-rj  sifiillnrly  ticrivetj  in  rofcrcncc  21.  It  is  notcw«rrfliy.  nion*ovfr.  that  snrh 
comliisittn.s  Imvc  also  bc»-n  cI»>riv(Hl  in  r»*f*‘rcnrcs  12  to  Hi  by  corisidfrnbly  dilTcrcnt  nii'thtKls  of 
nnnlysi.s.  'Phi*  {thysicfU  i‘xplnn'>,f ion  for  tbo  losscncil  flTcct  <if  fi  prt*s.siin*  j-Tiulimt  by  (smiirm 
of  the  wall  hns  hoi-n  ^'iven  In'lffH-nilcntly  in  rcfiTcncos  21  an<l  in. 


halii-a.  lo  daiiaasc  holli  lha  -k  m-l  ri)l  ii  ni  aoaHiaiiiil  .iml  lli. 
\us.st4l  numhar. 

.Silica  X  will  oriliiianlc  ha  mile  hilla  alfic  i iil  h\  ilu'  w.dl 
lamparalura.  acpiation  (21)1).  implio  llnil  loolna.:  of  ilia  wall 
will.  in  jraliaral,  land  lo  dimini>h  the  plic-iaal  hcinnd.iiv- 
later  ihic-kiii'ss  o.  llowiwer.  for  a  oivaii  t.-ihia  of  7  7  / 

tha  hoiindary-layar  ihiaknass  o  will  land  lo  nnicii'a  wii!. 
Mai-h  mimhar,  l•spaaially  in  a  favornhla  pra"mc  yradi.  ni 

sKi'AR  vrios 

For  a  fixed  valoaily  disl  rihni  ion  u  n  £-  onl-ida  ili, 
hoiindary  layer  and  a  lixad  .Maah  nninhar  .1/  .  dmiini'IniiL: 
till’ wall  ll■mpa|•al  ura  w  ill  land  lodi4a\  saparalioii  l)\  nioc  iiil: 

I  ha  sapara  I  ion  point  downsiraam.  'I’his  can  ha  ^aaii  parii.ii- 
larly  from  acpiation  (ti.'l).  aaaordina-  in  whiah  lha  \alna  of  ,\ 
rc'cpiircMl  for  saparalioii  X,,,,  w  ill  inaraa-.i'  a'  h  i-  cinnmi'had 
This  is  a  fiirthar  illiisi  rat  ion  of  tha  diminnlion  of  iha  cliiec  i 
I'fl'ai-I  of  a  prassiira  oradiaiil  liii  this  aasa.  an  adceisa  pii  '-nia 
eradicml)  hy  aoolinir  of  tha  wcill. 

'Phc'  idl'aal  of  .Mac4i  nnmhi'r  on  the  sc.paralion  ponii  for  a 
lixacl  dist  rihiil  ion  of  i/|  i/  (fl  and  cdl  hara  lixad  c  aliiaof/.  7  / 

ora  lixacl  value  of  l'„  1'  aaimol  hc‘ so  raadilt  predial  ad  from 
lha  acpialions  dl■valc)pacl  hara.  sinaa  an  inaraasa  of  .Maah 
numheu' in  l  hc>  advensa  prassurc' oradiaiil  w  ill  land  In  daaraasa 
hot  li  X„,,  (1X1,  (I).’! )  .1  and  X,,  ,,(4 1  ( aip  l()4  1 1 .  I  low  a\ er.  iiniiiariaal 
l•.xamplas  aarriad  out  for  lha  aasa  a,  a  I  J  ha\  c‘  indiaalad 
that  for  a  lixad  ralio  h  of  w  all  lainparal  lira  to  acpiilihrinm 
adiahatia  wall  lamparalura.  inaliidino  lha  aasa  of  an  insniataci 
wall  an  ini’raasc.  of  .Maah  iiiinihi'r  lands  to  anhanaa 

separation  hy  moviiifr  the'  saparation  point  npsiriaiin  laf. 
lalila  i  and  raf's.  2.  21,  III.  and  511).  I.owariiii;'  the  ili.xc'd 
valiia  of  //.  howc'var.  lands  to  diminish  this  imfavorahla 
all'i'al  of  Maah  numhar  on  saparation  (af.  raf.  21  i. 

If  lha  ratio  r„  '}\  of  wall  tampiu'al  ura  lo  fraa-st rc‘am  or 
rcd'ari-nac'  tamparaliira  inslc'iid  of  that  of  wall  tampi'rai  lira  to 
('(piilihrium  adiahatia  wall  ic'inparalura  h  is  ka))i  lixacl.  the 
afl'c'cl  of  Maah  numhar  on  tha  separation  ixiini  is  ahanirail 
'Phis  is  c'ssanlially  due  lo  tha  fac  t  that  for  a  lixacl  value  of 
/','7',  the'  lamparalura  ralio  h  claaraasas  with  Maah  nninhar 
(af.  i'(|.  (25) )  and  halic'i'  X,,,,  will  no  lonirc'r  tc'iid  to  ha  so  arc'al  l\ 
daarc'asaci  hy  an  inaraasa  in  (af.  c‘(|.  (loll.  ( 'onsacpiaiil  ly . 
tha  idl'cM’l  of  an  inarc'iisa  in  Maah  numhar  is,  in  iratiaral,  miiah 
lass  unfavorahia  in  this  aa.sa  and  ma,\ .  in  aartain  c  ase's,  move 
the-  saparation  point  dow  list raani.  aspaaiallx  at  hiydi  lixad 
values  of  7',,  7’,  . 

Fillin'  2  shows  the'  saparation  point  for  lha  aasa  a,  a 
1— t  as  a  function  of  lha  wall-1  amparal ura  ratio  7',,  7'  for 
.1/^  -  ().  'Pha  favorahlc'  afl'c'c-t  of  aoolinir  of  lha  wall  is 
alc'arly  si'c'ii  hara.  Fi<rura  4  shows  tha  si'])aral ion  point  as  a 
fimation  of  Maah  niimhc'r  for  a  lixc'd  ralio  of  wall  to  frc-c'- 
straam  rc'faraiiaa  Ic'inpi'ratura  T„  I\  2.  .Vn  inaraasa  of 
Maah  numhar  is  scan  in  this  aasa  aaliially  lo  move'  lha  sc'pa- 
ralion  ])oinl  downslrc'iim,  in  aontrasl  with  its  ali'aat.  also 
shown  in  lieiirc'  4,  at  zc'ro  lic'iit  transfc'r.  Fiirnras  2  and  4 
are  hased  on  t  ha  ('cpiations  dc'vc'lopad  hara,  and  furl  liar  dc'lails 
of  till'  aalatilations  aan  hc'  found  in  rafarc'iiaa  7. 

STABILITY  CHARACITERISTH'S 

It  has  already  haati  |)ointad  out  that  the  mat  hods  davalopad 
hara  may  hc'  c'xpc'atc'd  lo  yic'ld  siilliaiantly  aaaiirala  rasiills 
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I'lio  liK  2. — Srparalioii  point  as  a  function  of  wall  tcntpcratiirc.  In- 
coinprcssiLlc  flow:  k,  n  =  1  -  .U„  =  0. 


Fifn-RK  3. — Separation  point  as  a  function  of  Macli  inimber. 

—  1  -i. 

for  laminar-bouiularv-liiver  stability  calciilatioiis.  lii  fact, 
altliougli  stability  calctilatioiis  and  their  results  are  shown 
to  some  extent  in  referenees  1  and  2,  the  chief  purjxtse  of 
these  calculations  was  to  show  that  the  results  obtainetl  by 
the  approximate  methods  presented  here  compare  sufficiently 
closely  with  those  obtained  by  known  exact  solutions.  Such 
was  also,  at  first,  the  purpose  of  reference  .3.  Thus,  if  was 
found  that  for  compressible  flow  without  a  pressure  gradient, 
such  as  flow  over  a  flat  plate,  the  minimum  critical  Reynolds 
numbers  for  various  Mach  numbers  at  zero  heat  transfer, 
as  well  as  at  various  uniform  wall  temperature.s,  were  pre¬ 
dicted  with  satisfactory  accuracy  by  the  solutions  obtained 
by  the  methods  presented  here.  Moreover,  it  was  shown 


I  that  the  inaxiiuiim  wall  Icinpcrai  iii  i--  ito  l)i'  ciilliMl  hn  e  ihc 
“critical  temperatures'')  required  I,i  siablli/.e  llie  lluw  eoiu- 
I  pletely  were  also  ealeulaled  as  fuiielions  of  the  Mach  uum- 
j  ber  with  .satisfactory  tieeiiraey  on  the  l)asis  of  lliese  met  hods 
j  (presented  in  ttie  section  “Flow  Without  iin  .Vxial  l’re->utc 
dnnlwiit  With  .irbitrary  (( 'uiistiiiit  l  I’nniilll  .\ umber  an  ) 
]  Variable  Wall 'I'emperat  ttre"  I .  'I  he  itiitiitniitn  erit  ieal  K(  \ - 
Holds  number  for  ineoinprc'ssible  How  iti  tlu'  \  ieitiity  of  a  for¬ 
ward  stagnatiott  |)oint  its  caieuliiied  by  the  Kilt tmiu-i'oli!- 
haliscn  method  (cf.  the  siibseetion  ''.'s|  jigna t  ion  Flow  s"  in 
llic  scctioit  “F'low  ^Vdlh  Fresstire  ( iradleni ,  Praiidil  .\uiubei 
i  /’/•=- 1 .  and  rniform  Wall 'retiipeiitture  (A  (  otisiiint  i"  >  wn^. 
I  found  to  agree  well  with  ihiit  ealeulaled  liy  (heexiid  s<ihilion 
!  of  refenmee  a.'f.  Finally,  for  itieomjtressible  How  with  a 
i  linearly  tliminishitig  vidoeitx  outside  tlu>  hotiuilary  l;t\er.  the 
present  method  of  ciilculalion  (ef,  the  sei  lioii  “Flow  Wiib 
Fresstire  (iiiidietil.  Friiiidll  Siiiiilici'  /’r  I.  mill  l  iiiform 
Wall  'reinperattire  (A  t 'onstiitil  i"  i  w  its  found  to  lead  to  a 
minimitm  eritieal  Ri’ytiolds  number  in  sitt isfaet or\  agree- 
mc  It  with  that  ealciilated  from  tin'  solution  in  ri  feii  iiee  .",2 
I  Most  of  the  slabilil  \  calciilalious  which  hnvr  bi'eii  i-arried 
out  in  this  eoimtrv  have'  been  based  on  the  tituil\si-  iiiid 
criteria  developi'd  by  l,in  (ref.  .fS)  for  itieottipressible  How 
and  subse((uctltly  (>.\tetided  by  (at)  and  Lees  (refs  S  iitid  40 
to  cotn|)r<“,s.sible  How  ,  In  reference  S.  simplified  approximali' 
two-diimmsiotial  stability  eriteria  for  eompressihle  How  have 
been  develo|)etl.  whereby,  withotit  miieli  diHiiuilty.  it  is 
possible  to  caleulate.  for  a  givi'ti  type  of  How.  the  miiiimiim 
<-rilica)  Heynohls  numbers  as  well  tts  the  wttll  temperal  tiri' 
retptired  for  infiniti'  mitiimitm  l•riIieal  Ke>tiolds  number. 
As  will  be  c.\])lained  sttbsetpumlly,  these  criteria  lut\ i' 
recently  bi'eii  modified.  The  minimtim  eritieal  Ki'ynolds 
number  /f,,  „  is  the  minimum  Reynolds  ntimber  lu'eessary 
for  the  possibility  that  very  small  (list iirbanees  iti  the  boittid- 
ary  layer  may  be  amitlified  with  time;  that  is.  A',,  is  the 
minimum  Reynolds  ittmiber  riapiired  for  instability  of  the 
laminar  boundaty  layer  with  res|)eci  to  small  dist tirbaiiees 
of  at  h'ast  ccrtaiti  wavtdettgths.  The  wall  temperatitre  for 
infinite  values  of  AL  r,  is  tliett  itstiallv  ittterpreled  its  ilie 
highest  Icmperattue  foi'  which  the  (compressible)  laminar 
boundary  layer  will  be  cmnplett'ly  stable  for  all  Ri'vtiolds 
numbers.  The  analyses  in  refcretiees  .fS.  40.  and  S  atid 
subseciuent  analyses  based  on  them  are  of  practical  interest, 
since  under  the  eonditioti  of  a  suflieietttly  low  free-streain 
turbiih'iicc  a  nece.s.sarv  (though  not  stiflieient)  eotidition  for 
transition  from  a  lamitiar  to  a  turbtilettl  botitidary  layer 
appears  to  be  an  instability  of  the  laminar  la  wr.  A  siirviy 
(as  of  19.52)  of  tlu'oretieal  and  cxperinnuital  investigations 
on  laminar-boundary-layer  stability  can  be  found  in  reference 
.5!). 

Tile  purpose  of  tlie  present  sub.seclion  is  to  summarize  the 
theoretical  investigations onlaminar-boundarv-layer  stability 
performed  at  the  Folytechnic  Institute  of  Brooklyn  by  using 
the  mean-flow  (or  steady-state)  solutions  obtained  by  the 
mctfiods  prc.senicd  in  this  n'port.  In  an  unpublished  ri'port 
entitled  “( 'alculation  of  Stability  of  ( 'onstant-Fressure 
Boundary  Layers  on  Isothermal  Surfac('s  With  an  Integral- 
Method  Mean-Flow  Solution”  Frofessor  Martin  Bloom  dc- 
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vi'lopi'd  (■(‘rtaiii  nioditiculioiis  of  Lees’  approximate  stability 
eriteria  (ref.  8)  and  applied  these  to  the  ealeulation  of  the 
stal)ility  of  the  laminar  homidary  layer  over  a  flat  plate  at 
uniform  wall  temperature.  'Phis  report  is  available  for  loan 
or  reference-  in  the  JOivision  of  Research  Information,  Na¬ 
tional  Advisory  (’ommittee  for  Aeronautics,  Washington, 

1) .  (,’.  'Phis  work  is  summarized  in  references  -i  to  5. 

Minimum  critical  Reynolds  numbers  for  given  wall  tempera¬ 
tures  and  Mach  numbei-s  were  calculated.  Moreover,  the 
wall  temperature  required  to  stabilize  the  flow  completely 
was  also  calculated  as  a  function  of  the  Mach  number. 
Similar  types  of  calculations  with  similar  results  were  carried 
out  independently  by  Van  Driest  (ref,  60),  and  these  are  now 
well  known."  Briefly,  the  results  indicate  the  stabilizing 
effect  of  cooling  of  the  wall  by  increasing  the  minimum 
critical  Reynolds  number  for  a  given  Mach  number.  More¬ 
over,  for  a  Prandtl  number  Pr  of  1,  it  is  foimd  that  the 
boundary  layer  can  be  completely  stabilized  by  sufficiently 
low  wall-temperature  ratios  for  .Mach  numbers  .!/„ 

between  1  and  approximately  o.  (For  /V=().72,  this  can  be 
tlieoretically  accomplished  for  1  At  higher 

Mach  numbers,  particularly  in  the  hypei-sotuc  range,  the 
validity  of  the  theoretical  approach  has  not  been  establisbetl. 

The  stability  of  the  laminar  compressible  boundary  layer 
in  a  pressure  gradient  has  been  analyzed  iti  n-ference  44  for 
zero  heat  transfer  at  the  wall.  Calculations  there  for  the 
supersonic  flow  over  a  thin  biconvex  airfoil  indicated  the 
stabilizing  influence  of  the  favorable  pre.ssure  gradient. 
Tliis  stabilizitig  influence,  however,  was  found  to  be  con¬ 
siderably  diminished  at  higher  free-stream  .Mach  numbers 
M^=A.  The  stabilizing  influence  of  a  favorable  pressure 
gradient  cati  also  be  clearly  illustrated  by  comparing  the 
minimum  critical  Reynolds  number  (namely, 

2.40X10”  (ref.  2))  for  the  incompressible  flow  =  f  in 
the  vicinity  of  a  forward  stagnation  point  with  the  much 
smaller  value  /i’a,.cr|=7.3X  10’  (ref.  1)  for  incompressible 
flow  over  a  flat  plate.  The  destabilizing  effect,  in  the  case 
of  zero  heat  transfer,  of  an  adverse  pressure  gradient  is 
readily  illustrated  by  considering  the  case  uju^  =  l—^  (ref. 

2) .  The  minimum  criti<-al  Reytiolds  numbers  for  this  case 
for  .1/^=0  and  1  are  compared,  in  table  II,  with  the  larger 
values  for  flow  over  a  flat  plate  taken  from  reference  1 . 

For  compressible  flows  with  heat  transfer  and  pressure 
gradient,  the  only  stability  calculations  which  appear  to  have 
been  made  thus  far  are  those  in  references  7  and  62.  In 
both  of  these  references  oidy  Mach  numbers  of  3  or  lower 
were  considered.  (C’f.  footnote  11.)  In  reference  62,  the 
small-perturbation  solutions  of  reference  25  are  used,  while 
reference  7  uses  solutions  based  on  the  methods  of  analysis 
presented  in  the  present  report.  Reference  62  shows  that 

’•  Hloom’.s  calciilaiions  (refs.  anti  4;  jjaw  r'*sults  quite  similar  l«  the  wpll-kiiown  n»- 
suUs  t>f  Van  Driest  (ref.  ffO).  Further  modlfifatlons  of  the  stability  crlterih.  however,  Ie«l  to 
rather  ef>mpllrart‘fl  eurve.s  with  several  brariehes  of  critical  temperature  raUt>  versus  Mach 
number  (ref.  Thes«>  were  apparently  due  tt)  larce  values  of  the  stability  parameter  X  'not 
reMfetl  to  the  X  of  the  piesent  rpjiort)  as  deflntMl  in  referenei*  s.for  tarve  va!ue.s(>f  Af*  or  small 
values  of  7»/7’(5e.  Dunn  and  Lin  (ref.  (iIi,  howt»ver,  have  tjuite  recently  made  basic  refine¬ 
ments  in  the  analysis  of  reference  40  anti  have  develop<*d  a  more  accurate  set  of  both  two- 
dimensional  anti  three-tlinienslonal  stability  criteria.  Calculations  for  flow  over  a  Hat  plate 
hasetl  on  this  set  yieldetl  results  qtilte  similar  to  those  of  Van  Driest  orof  1llf>om*sflrsf  ealcula- 
tloas.  Aecordinp  to  the  Diinn-Lin  erileiia,  the  values  of  X  (now  retlefinetl)  remaimHl  quite 
small  for  flow  over  a  flat  plate  even  at  hl!?h  Mach  numbers.  The  new  two  iIimeDslunal 
criteria  do  not  ap|H*ar  to  yield  results  aprireclably  tllfferent  from  tho.se  of  Bloom  or  Uvs  for 
Mach  numbers  below  appnixlmately  :i. 


Fidj  rk  4. — Uarios  of  wall  to  local  inaiu-siroani  ictiiixTatiirc  rotiuirocl 
for  Infinite  ininiinnni  critical  Kcynold.s  niiniluT  ov(‘r  thin  biconvi'x 
supersonic  airfoil. 


the  critical  wall  temperatures  rcepiircd  to  stabilize  the  lami¬ 
nar  boundary  layer  completely  arc,  for  a  given  Mach  num¬ 
ber,  increased  by  a  favorable  pressure  gradient  and  decreased 
by  an  adverse  pressure  gradient.  Further  calculations  also 
indicate  the  greater  amount  of  cooling  required  to  stabilize 
completely  the  flows  with  adverse  pressure  gradients  than 
that  required  for  those  with  favorable  pressure  gradients. 
This  illustrates  in  a  further  fashion  the  stabilizing  influence 
of  a  favorable,  and  the  destabilizing  influence  of  an  adverse, 
pressure  gradient.  In  reference  7,  critical  wall  temperatures 
have  been  determined  for  the  supersonic  flow  over  a  thin 
biconvex  airfoil  at  two  given  stations  along  tlic  flow,  and 
these  have  been  compared  with  the  corresponding  results 
for  flow  over  a  flat  plate.  The  results  are  shown  in  figure  4. 
wherein  it  is  seen  that  higher  critical  temperature  ratios 
T„lTi  are  obtained  for  the  flow  with  the  favorable  pressure 
gradient  than  for  the  flow  over  a  flat  plate.  It  may  be  ob¬ 
served,  in  this  eonneetion,  that  for  a  given  reference  tem¬ 
perature  Te,  at  a  point  immediately  behind  tlic  shock  wave 
at  the  leading  edge  of  the  supersonic  airfoil,  the  critical  wall 
temperature  may,  at  the  higher  Mach  numbers,  be  greater 
for  the  favorable-pressure-gradient  case  than  for  the  flat- 
plate  ease.  This  is  due  simply  to  the  fact  that  the  local 
temperature  Tt  outside  the  boundary  layer  over  the  airfoil 
diminishes  along  the  flow  (see  ref.  7  for  details). 


ANALYSIS  OF  LAMINAR  COMPRESSIBLK  BOrNDARY'  LAYER  WITH  ilEAT  TRANSFER  AM)  PHFSSI  RF  (iHADIF.NT 


17 


FKii  HK  5. — Miniiimiu  (Tifical  Reynolds  immLcrs  Wt,  c,  vcrsii.v  Mai-h 
iiutnLcr  .U4,  wall-rcforeiicp-t(>nip<Tature  ratio  T„/Ti,  and  wall- 
ftiiiilibrium-tpinpcratiirp  ratio  h  Vortical  lines  with  arrows  indi¬ 
cate  corresponding  asymptotes;  compare  with  tattle  III. 

A  second  type  of  stability  calculation  carried  out  in  ref¬ 
erence  7  was  the  determination  of  minimum  critical  Rey¬ 
nolds  numbers  for  the  laminar  boundary  layer  at  a  given 
station  of  the  supersonic  airfoil  for  various  values  of  the 
(uniform)  wall  temperature  and  free-stream  Mach  number. 
The  results  are  shown  in  table  III  and  figure  5,  where  com¬ 
parison  is  also  made  with  the  flow  over  a  flat  plate.  The 
stabilizing  effect  of  cooling  of  the  wall  and  of  the  favorable 
(negative)  pressure  gradient  here  can  be  clearly  seen. 

From  table  III  and  figure  5,  the  effect  of  Mach  number  on 
the  stability  characteristics  is  seen  to  depend  on  the  pressure 
gradient  and  whether  the  ratio  A=7'„/7’t  of  wall  tempera¬ 
ture  to  equilibrium  adiabatic  wall  temperature  or  the  ratio 
TJTt,  of  wall  temperature  to  reference  temperature  is  held 
fi.xed.  PTom  figure  5  it  is  seen  that  for  a  fixed  h  an  increase 
of  Mach  number  from  1.5  to  2.0  destabilizes  the  boundary 
layer  both  over  a  flat  plate  and  over  the  airfoil.  This  effect 
is  seen,  in  fact,  to  be  enhaiu-ed  by  the  negative  pressure 
gradient  here.  For  a  fixed  value  of  the  ratio  TJ1\,  how¬ 
ever,  an  increase  of  Mach  number  is  now  seen,  from  figure 
5,  to  have  a  stabilizing  influence  on  the  flow  without  a 
pressure  gradient,  especially  at  the  lower  wall  tempera¬ 
tures.  For  the  flow  over  the  airfoil,  however,  figure  5  (cf. 
also  table  111(a))  now  indicates  that  an  increase  of  Mach 
number  has  a  stabilizing  effect  only  at  wall  temperatures 


I  close  to  the  critical  temperature  and  that  for  (fixed  1  liiofier 
I  wuli-temja'rat lire  ratios  of  7',,  /),  an  iucnaise  of  Mach  num- 
j  her  has  a  clear  di'staliilizing  effect  similar  to  the  ease  (d 
fixed  //. 

(  ()N<  LL.SI()N.S 

From  the  analysis  of  eompressihle  laminar  hoimdarx 
layers  with  lieat  transfer  and  with  and  without  pic'siire 
gradient  pri'senttal  herein  ittider  the  assntnptioti  of  a  linear 
■  temperature-viscosity  relation,  the  following  comhisiori' 
can  be  ilrawn : 

1.  P'or  flow  without  a  prissiire  gradient,  such  as  flow  o\ei 
a  flat  plate,  the  botindary-layer  eharaeteristies  can  he  easih 
delermineil  from  the  eipiatiotis  developed  here  for  a  given 
constant  I’randtl  tuimber  (of  the  order  of  iinigtdinde  id' 
j  unity  ),  a  given  Mach  iitimber,  and  a  given  «  all-temperat  ore 
I  distribution. 

,  2.  For  flow  with  a  pressure  gradient,  the  houndary-layei 

I  characteristics  can  also  be  easily  determined  from  the 
I  erpiations  developed  here,  provided  the  I’raiidtl  miinher  i^ 
I  unity  and  the  wall  temiieratiire  is  uniform.  Here,  the  \  e- 
locity  distribtitiot)  otitside  the  botitnlary  layer  atid  the  frei'- 
slream  Mach  number,  as  well  as  the  wall  temperature,  are 
considered  as  prescribed.  'I'he  etpiatiotis  are  also  valid  for 
zero  h(>al  tratisfer  at  the  wall  {h  -  I  where  h  is  the  ratio  of 
stagtialion  etithalpy  at  the  wall  to  stagnation  etithal|)y  at 
•he  outer  edge  of  the  botitnlary  layer). 

d.  A  ndatively  sim|)lc  method  of  calciilat  ing  the  separat  ion 
point  iti  a  givtui  stibsonic  or  sttpersotiic  ad  terse  pressure 
gradient  over  a  wall  at  any  spi'citied  tmiform  lemperatttre 
lias  b<*<'ii  dcvt'Joped  Iwiv.  'I'his  imdhoil  is  also  applicable  for 
zero  heat  transfer  {//—■  I ). 

4.  A  comparison  of  the  results  of  the  iindhods  in  coti- 
clusions  1,  2.  and  d  with  known  exact  sohitions  for  various 
types  of  flows  indicates  that  the  nndhods  of  ealctilation 
developed  here  may  be  expt'cted.  in  general,  to  yield  rt'stilts 
of  sufficient  aoctiracy  for  |)raetieal  ptirposes.  incltiding  sta¬ 
bility  calculations. 

5.  P’rom  the  equations  developed  here,  it  can  In-  shown 
that  coolitig  of  the  wall  tetnls  todimitiish  the  .Xtisselt  ntimber 
and  especially  the  skiti-frictioti  coefiieii'tit  iti  a  favorable 
(negative)  pre.ssure  gradient  and  to  iticrease  the  coefficietits 
in  an  adverse  pressure  gradietit.  Hecattse  of  the  |)ropor- 
tionality  factor  iti  the  viscosity-temiierature  relation  assumed 
here,  it  also  follows  that  lowering  the  ratio  of  wall  to  free- 
stream  temperature  will,  inde])end('ntly  of  the  pressure 
gradient,  ordinarily  tend  to  im'reasc  both  the  Nussidt  number 
aiul  the  skin-friction  coefficient. 

6.  The  equations  developed  here  further  inqily  that  cool¬ 
ing  of  the  wall  tends,  in  general,  to  diminish  th(>  direct  effect 
of  a  pressure  gradient,  while  heating  tends  to  enhance  it. 
A  particularly  clear  example  of  this  is  the  delay  of  si'paration 
in  an  adverse  pressure  gradient  by  cooling  of  the  wall. 

7.  The  results  of  a  numerical  example  for  a  fixed  linearly 
decreasing  velocity  outside  the  boundary  layer  indicate,  in 
addition  to  the  delaying  of  separation  by  cooling  of  the  wall, 
that  for  a  fixed  ratio  h  of  wall  temperature  to  e((uilibrium 
adiabatic  wall  temperature  an  increase  of  free-stream  Mach 
number'moves  the  separation  point  upstream,  while  for  a 
fixed  ratio  of  wall  temperature  to  free-stream  ttunperature 


IS 
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t'a'  R'l  iiUTeasc  of  Mach  iiumhcr  has,  in  jiciuTal,  a  less 
imfavorahlc  cfFcct  and  in  tlic  ease  7’,, '7'^  =  2  actually  moves 
the  seiianition  point  downstream. 

S.  M'hih“  cooling  of  the  wall  tends,  in  j;eneral,  lo  stabilize 
the  laminar  houndarv  layer,  it  is  shown  theoretically  that  at 
in  xlerate  supersonic  Mach  iiiimhers  sufficient  cooling  may 
completely  stabilize  the  boundary  layer.  At  higher  Maidi 
numbei's,  particularly  in  tlie  hypersonic  range,  the  validity 
of  the  theoretical  approach  has  not  been  established.  A 
favorable  pressure  gradient  has,  in  gi'iicral,  a  stabilizing 
effect  on  the  laminar  boundary  layer,  while  an  advei’se 
pressure  gradient  has  a  destabilizing  effect.  A  numerical 
e.xample  for  supersonic  flow  over  a  thin  airfoil  illustrates  in 
detail  these  and  other  effects  of  .Mach  number,  wall  tem- 
|)erature,  and  pressure  gradic'iit  on  the  stability  of  the  ! 
laminar  boundary  layi'r.  I 
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